Dynamical Dark Matter: I. Theoretical Overview 
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In this paper, we propose a new framework for dark-matter physics. Rather than focus on one or 
more stable dark-matter particles, we instead consider a multi-component framework in which the 
dark matter of the universe comprises a vast ensemble of interacting fields with a variety of different 
masses, mixings, and abundances. Moreover, rather than impose stability for each field individually, 
we ensure the phenomenological viability of such a scenario by requiring that those states with larger 
masses and Standard-Model decay widths have correspondingly smaller relic abundances, and vice 
versa. In other words, dark-matter stability is not an absolute requirement in such a framework, 
but is balanced against abundance. This leads to a highly dynamical scenario in which cosmological 
quantities such as JIcdm experience non-trivial time-dependences beyond those associated with the 
expansion of the universe. Although it may seem difficult to arrange an ensemble of states which 
have the required decay widths and relic abundances, we present one particular example in which 
this balancing act occurs naturally: an infinite tower of Kaluza-Klein (KK) states living in the bulk 
of large extra spacetime dimensions. Remarkably, this remains true even if the stability of the KK 
tower itself is entirely unprotected. Thus theories with large extra dimensions — and by extension, 
certain limits of string theory — naturally give rise to dynamical dark matter. Such scenarios also 
generically give rise to a rich set of collider and astrophysical phenomena which transcend those 
usually associated with dark matter. 



INTRODUCTION, MOTIVATION, AND 
SUMMARY 



Situated at the nexus of particle physics, astrophysics, 
and cosmology lies one of the most compelling mysteries 
that faces physics today: that of unravelling the identity 
and properties of dark matter From measurements 
of galactic rotation curves and velocity dispersions to ob- 
servations of the gravitational lensing of galaxy clusters 
and the detection of specific acoustic peaks of the cosmic 
microwave background (CMB), ample circumstantial ev- 
idence suggests that most of the matter in the universe 
does not interact strongly or electromagnetically. Such 
matter is therefore electrically neutral (dark) and pre- 
sumed non-relativistic (cold). Beyond these properties, 
however, very little is known about the nature of dark 
matter. Fortunately, the current generation of dark- 
matter experiments have unparalleled sensitivities, and 
new data concerning the possible direct and indirect de- 
tection of dark matter can be expected soon. This data 
will therefore go a long way towards not only resolving 
this pressing cosmological mystery, but also constraining 
the possibilities for physics beyond the Standard Model 
(SM). 

Many theoretical proposals for physics beyond the 
Standard Model give rise to suitable dark-matter can- 
didates. However, most of these dark-matter candidates 
consist of a single particle (or a small collection of par- 
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tides) which are stable on cosmological time scales as 
the result of a discrete symmetry. Examples include 
the lightest supersymmetric particle (LSP) in supersym- 
metric theories, and the lightest Kaluza-Klein particle 
(LKP) in certain higher-dimensional theories in which 
the Standard Model propagates in the bulk [2|. In the 
first case, the LSP is stabilized by the assumption of an 
i?-parity symmetry, while in the second case the stabi- 
lizing symmetry is a so-called "KK parity" . However, 
in all cases, the ability of these particles to serve as 
dark-matter candidates rests squarely on their stability. 
Indeed, any particle which decays into Standard-Model 
states too rapidly is likely to upset traditional big-bang 
nucleosynthesis (BBN) and its successful predictions of 
light-element abundances. Such decays can also leave un- 
acceptable imprints in the cosmic microwave background 
and diffuse X-ray and gamma-ray backgrounds. For this 
reason, stability is often the very first criterion required 
for the phenomenological success of a hypothetical dark- 
matter candidate. 

There is, of course, one important exception to this ar- 
gument: a given dark-matter particle need not be stable 
if its abundance at the time of its decay is sufficiently 
small. A sufficiently small abundance ensures that the 
disruptive effects of the decay of such a particle will be 
minimal, and that all constraints from BBN, the CMB, 
etc., will continue to be satisfied. 

In this paper, we shall consider a new framework for 
dark-matter physics which takes advantage of this possi- 
bility. Specifically, we shall consider a multi-component 
framework in which the dark matter of the universe com- 
prises a vast ensemble of interacting fields with a variety 
of different masses, mixings, and abundances. Rather 



2 



than impose stability for each field individually (or even 
for the ensemble of fields as a whole) , we shall ensure the 
phenomenological viability of such a scenario by requir- 
ing that those states with larger masses and larger decay 
widths into Standard-Model fields have correspondingly 
smaller relic abundances, and vice versa. In other words, 
dark-matter stability is not an absolute requirement in 
such a framework, but is balanced against abundance. 
As we shall demonstrate, this leads to a highly dynam- 
ical scenario in which cosmological quantities such as 
f^CDM experience non-trivial time-dependences beyond 
those associated with the expansion of the universe. We 
shall therefore refer to such a scenario as "dynamical dark 
matter" . 

In general, it might seem difficult (or at best fine- 
tuned) to have an ensemble of states which are not only 
suitable candidates for dark matter but in which the 
abundances and decay widths are precisely balanced in 
this manner. However, it turns out that theories with 
large extra spacetime dimensions not only naturally pro- 
vide such ensembles of states, but do so in a manner 
which is virtually intrinsic to their construction. If the 
Standard Model is restricted to a brane floating in a 
higher-dimensional space, it then immediately follows 
that any field propagating in the bulk of this space must 
be neutral under all Standard-Model symmetries. As a 
consequence, such bulk fields can have at most gravita- 
tional interactions with the physics on the brane, and will 
therefore appear as dark matter from the perspective of 
an observer on the brane. Moreover, from the perspective 
of this four-dimensional observer, such bulk fields will ap- 
pear as an infinite tower of individual Kaluza-Klein (KK) 
modes. This, then, would constitute our dark-matter en- 
semble. 

At first glance, such a scenario for dark matter would 
appear to face a major phenomenological hurdle: in the 
absence of additional symmetries or ad-hoc assumptions, 
an entire Kaluza-Klein tower of bulk states will generally 
be unstable: the heavy KK states in such a tower will 
generically decay into not only lighter KK bulk states 
but also Standard-Model brane states, and the lighter 
KK states will also decay into Standard-Model brane 
states. Even the stability of the lightest modes of the bulk 
field is not guaranteed. This instability of the Kaluza- 
Klein tower therefore appears to pose a serious threat for 
the survival of big-bang nucleosynthesis in its traditional 
form, and can similarly disturb the X-ray and gamma-ray 
backgrounds. 

Fortunately, there are two critical features of Kaluza- 
Klein towers which can play off against each other in or- 
der to render such a scenario phenomenologically viable. 
As one goes higher and higher in a generic KK tower, it 
is true that the decay width of the KK states generally 
increases with the KK mass. However, it is also true that 
the cosmological abundance associated with such states 
can often decrease with the KK mass. This is particularly 
true if we imagine that these states are cosmologically 
produced through misalignment production, as turns out 



to be particularly appropriate for such scenarios. As a re- 
sult, it might be possible that all KK states which decay 
before or during BBN have such small abundances that 
the destructive effects of their decays are insignificant, 
while at the same time a significant fraction of the KK 
tower survives to the present day and thereby contributes 
to the observed total dark-matter abundance. Thus, the 
surviving dark matter at the present day would consist 
of not merely one or two states, but a significant fraction 
of an entire interacting KK tower. Through the existence 
of such "dark towers" , theories of large extra spacetime 
dimensions therefore provide an ideal realization of our 
general dynamical dark-matter scenario. 

In this paper, we shall lay out the general properties 
of such a scenario and explore the extent to which such 
a scenario is viable. Moreover, we shall attempt to do so 
in a completely model-independent way, without mak- 
ing any assumption concerning the nature of the bulk 
field in question. However, it is important to recognize 
that this entire approach represents a somewhat unortho- 
dox approach to dark-matter physics. By balancing the 
stability of the different dark-matter components against 
their abundances across a large or even infinite ensemble, 
the dark matter in this scenario is intrinsically dynamical 

- its different components continue to experience non- 
trivial mixings and decays throughout their cosmological 
evolution, with such dynamical behavior continuing un- 
til, during, and beyond the current epoch. Moreover, be- 
cause the dark matter in our scenario has multiple com- 
ponents, its phenomenology cannot be characterized in 
terms of a single mass or annihilation cross-section. This 
can therefore lead to an entirely new dark-matter phe- 
nomenology, profoundly changing the way in which such 
dark matter might be observed and constrained through 
collider experiments and astrophysical observations. In- 
deed, within the specific context of large extra dimen- 
sions, we shall see that one important new phenomenon 
that emerges for such dark matter is the possibility of 
"decoherence" — i.e., the phenomenon in which only a 
single linear combination of KK modes couples to brane 
physics at one instant before decohering and becoming 
essentially invisible to the brane at all subsequent times. 

It is also important that this framework not be con- 
fused with recent proposals concerning so-called "Kaluza- 
Klein dark matter" Q ■ Theories of KK dark matter re- 
quire that the entire Standard Model propagate in the 
large extra dimensions and that the lowest ex- 

cited KK mode of a Standard-Model field (such as the 
lowest-excited KK photon or neutrino) be stable as the 
result of an internal geometric symmetry such as KK 
parity 0, H[ . Such theories of KK dark matter are there- 
fore similar to theories of supersymmetric dark matter 

— they are theories of a single, stable, dark-matter par- 
ticle. While phenomenologically consistent, such a point 
of view is diametrically opposed to what we are suggest- 
ing here. Moreover, in doing away with the infinite tower 
of KK states and focusing exclusively on the single light- 
est KK mode, such "KK dark matter" theories also do 
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away with that part of the physics which is intrinsically 
higher-dimensional. The resulting scenarios are there- 
fore insensitive to the rich physics that can emerge from 
an entire tower of Kaluza-Klein states acting in unison, 
with non-trivial masses and mixings governing their dy- 
namics. Indeed, it is precisely such behavior that would 
give a window into the nature of the extra dimensions 
from which such states emerge. 

By contrast, because our scenario balances a spectrum 
of decay lifetimes against a spectrum of relic abundances, 
our framework is sensitive to the physics of the entire 
tower of KK states and thereby makes use of the full 
higher-dimensional "bulk" in a fully higher-dimensional 
way. Moreover, since Type I string theories naturally give 
rise to closed-string states (such as the graviton, various 
moduli, and axions) which live in more spacetime dimen- 
sions than the Standard-Model open-string states which 
are restricted to live on D-branes, the scenario we shall 
be investigating is also extremely natural — and indeed 
almost unavoidable — in certain limits of string theory. 
Our work can therefore be seen as providing a test of 
the extent to which such string theories remain cosmo- 
logically viable as a function of the volume of the extra 
dimensions transverse to the Standard-Model brane. In 
other words, by studying dynamical dark matter and its 
phenomenological viability, we are not only exploring a 
new candidate for dark matter but also providing new 
phenomenological constraints on large extra dimensions 
and certain limits of string theory. 

This paper is organized as follows. In Sect. II, we shall 
introduce our dynamical dark-matter framework in its 
most general form, without making reference to the spe- 
cific example of large extra dimensions or KK towers. We 
shall discuss how lifetimes and abundances can play off 
against each other in such a scenario, and sketch the re- 
sulting contributions to the total dark-matter abundance 
as they evolve in time. In Sect. Ill, we shall then focus 
on the example of a generic tower of Kaluza-Klein states 
emerging from the bulk of large extra dimensions, and 
show that such a tower has all the required properties to 
be a dynamical dark-matter candidate, with abundances 
and lifetimes that satisfy unique mathematical inverse 
relations. We shall then proceed, in Sect. IV, to discuss 
several laboratory and astrophysical signatures of such 
a scenario, focusing on those new features which tran- 
scend typical dark-matter signatures and which might 
explain why such dark matter has not yet been observed. 
Throughout this paper, our analysis shall be as general as 
possible without specifying the precise nature of the bulk 
field in question. We shall then conclude in Sect. V with 
a discussion of extensions and possible generalizations of 
our dark-matter framework. 

This paper is the first in a two-part series. The primary 
purpose of the present paper is merely to provide a gen- 
eral theoretical overview of the dynamical dark-matter 
framework. As a result, we will not choose a particular 
species of dark-matter field, neither restricting ourselves 
to specific numbers nor subjecting ourselves to specific 



phenomenological bounds. Instead, our discussion here 
will focus on the full range of theoretical possibilities af- 
forded by this new scenario. However, in a companion 
paper Q we will provide a detailed "proof of concept" 
by focusing on the particular case that the KK bulk field 
in question is an axion. In particular, in Ref. [6| we will 
demonstrate that a bulk axion field can satisfy all the- 
oretical and numerical constraints needed to serve as a 
dynamical dark-matter candidate, and moreover we will 
demonstrate that this candidate also satisfies all known 
cosmological, astrophysical, and collider bounds on dark 
matter. Indeed, in making this assertion, Ref. [6j will 
borrow heavily from the results of a third paper 7] , a de- 
tailed forthcoming phenomenological study of axions and 
axion-like particles in higher dimensions. Thus, taken to- 
gether, these papers will demonstrate that our dynamical 
dark-matter scenario remains a very real possibility for 
explaining the dark matter of the universe. 



II. DYNAMICAL DARK MATTER: GENERAL 
SCENARIO 

In this section we shall begin by discussing our dy- 
namical dark-matter scenario in its most general form, 
without reference to specific examples such as those in- 
volving extra spacetime dimensions or KK towers of bulk 
fields. 



A. r versus Q,: Balancing lifetimes against 
abundances 

Broadly speaking, upon positing any new scenario for 
dark matter, one faces certain immediate constraints 
which must be satisfied. These constraints ultimately 
restrict either the abundance of the dark matter, the life- 
time of the dark matter, or the relation between the two. 

Let us begin by considering the case of a single dark- 
matter particle x- Since this dark-matter particle is 
presumed unique, it alone must carry the entire ob- 
served dark-matter abundance: il x — JIcdm ~ 0.23 0. 
However, given this large abundance, consistency with 
constraints coming from big-bang nucleosynthesis, the 
cosmic microwave background, and diffuse X-ray and 
gamma-ray backgrounds together require that \ have a 
lifetime which meets or exceeds the current age of the 
universe. Otherwise, decays of x run the risk of disturb- 
ing BBN and its successful predictions for light-element 
abundances. Such early decays also have the potential to 
distort the cosmic microwave background as well as the 
X-ray and gamma-ray backgrounds. 

However, because of the quantum-mechanical nature 
of the decay process, not all dark matter will decay at 
once. As a result, the lifetime of x must actually exceed 
the age of the universe by at least one or two orders 
of magnitude in order to ensure that x has a negligible 
chance of having already decayed in the recent past. This 
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likewise implies that such a particle also has a negligible 
chance of decaying either today or tomorrow. Such a 
particle x is therefore "hyper-stable" . Indeed, this is the 
case for most if not all known single-particle dark-matter 
candidates. 

Hyper-stability is the only way in which a single- 
particle dark-matter candidate can satisfy the compet- 
ing constraints of having a significant abundance Cl x <~ 
^cdm while simultaneously avoiding the dangerous ef- 
fects of decaying into Standard-Model particles. This 
then results in a dark-matter scenario which is "frozen" 
in time, with cosmological quantities such as fl x evolving 
only because of the expansion of the universe. 

However, the primary purpose of this paper is to pro- 
pose that there is another way — a "dynamical" way — 
to satisfy these competing constraints. First, we recog- 
nize that in some sense, it is somewhat unnatural to con- 
sider the dark matter of the universe as consisting of only 
a single particle. After all, the visible matter of the uni- 
verse constitutes only a small fraction of the energy den- 
sity attributed to dark matter, and yet is teeming with a 
diversity and complexity, known as the Standard Model, 
in which a complex network of elementary particles is 
organized according to their own internal principles. It 
therefore seems natural to consider the new opportunities 
that are open to us by taking the dark matter to consist 
of multiple particles as well. 

In this paper, we shall therefore imagine that the 
dark matter consists of a vast ensemble of particles \i-> 
i = 1, 2, N with N > 1. The first observation that fol- 
lows from this assumption is that none of these particles 
individually needs to have a significant abundance, since 
they may still collectively yield the correct total abun- 
dance ^cdm- As a special case, for example, we might 
imagine that each particle \i shares a common abun- 
dance fli = Hcdm/N. However, if these particles were 
also to have equal lifetimes, then this would not solve our 
second constraint — that of protecting the successes of 
BBN and minimally impacting the CMB and other dif- 
fuse backgrounds — unless each of these particles is not 
only stable but hyper-stable. This is because the net ef- 
fect of the nearly simultaneous decays of each of these N 
particles would be no different than that of the decay of 
a single particle carrying the full abundance ficDM- 

However, the fact that we have multiple particles fur- 
nishes us with an alternate way to satisfy these con- 
straints: we can imagine that each of these particles has a 
significantly different lifetime. In general, these particles 
can also have different individual abundances. As long 
as those particles which have relatively short lifetimes 
also have correspondingly small abundances, and as long 
as those particles which have relatively large abundances 
also have relatively long lifetimes, we can reproduce the 
correct total dark-matter abundance ficDM while simul- 
taneously avoiding any damaging effects on BBN, the 
CMB, etc. In this way, the existence of a vast ensem- 
ble of dark-matter particles \i opens up the possibility 
of balancing abundances against decay widths in a non- 



trivial way across a multitude of states. 

This, then, is the essence of our dynamical dark-matter 
proposal. The fact that we have distributed the to- 
tal required dark-matter abundance across many states 
means that no particular state is forced to carry a sig- 
nificant abundance on its own. We thus have the room 
to give these states a whole spectrum of lifetimes (or de- 
cay widths) without running afoul of cosmological con- 
straints. 

Note that the usual scenario of a single hyper-stable 
dark-matter particle is nothing but a special case of this 
more general framework: even though our scenario has 
N 3> 1, it still remains possible that one particle (or just 
a few particles) could carry the bulk of the abundance 
Ocdm at the present time. Following the same logic 
that applied in the single-particle case, this small sub- 
set of states would then be required to be hyper-stable, 
and all of the remaining states would have abundances 
that are far too insignificant to be of consequence. How- 
ever, the novel features of our scenario emerge in the 
opposite limit, when we imagine that none of our dark- 
matter states individually carries the bulk of the total 
relic abundance. Some fraction of these states then no 
longer need to be hyper-stable, leading to a dynamical 
scenario in which spontaneous dark-matter decays are oc- 
curring prior to, during, and beyond the current epoch. 
As a result of this behavior, cosmological quantities such 
as Ocdm will experience time- variations which transcend 
those due to the ordinary expansion of the universe. 



B. Dynamical dark matter: Time-evolution of 
individual components 

It is possible to outline the salient features of this sce- 
nario somewhat more quantitatively without loss of gen- 
erality. For this purpose, we shall describe the history 
of the universe as progressing through four distinct eras 
respectively associated with inflation (vacuum domina- 
tion), reheating (RH), radiation domination (RD), and 
matter domination (MD). Note that the reheating era is 
itself essentially matter-dominated, with the matter in 
this case consisting of the oscillating inflaton. Likewise, 
note that even though the current epoch is technically 
a ACDM universe, approximating this epoch as matter- 
dominated is nevertheless a fairly good approximation. 
For the purposes of our general discussion, we shall not 
specify the particular energy or temperature scales as- 
sociated with the transitions between these eras; such 
scales, especially as they relate to inflation and reheat- 
ing, are likely to be highly model-dependent. However, 
regardless of these energy scales, a quantity whose energy 
density p and pressure p are related through a single- 
parameter equation of state of the form p — wp will have 
a relative relic abundance f2 = p/p C mt that scales as a 
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function of cosmological time t according to 1 

f i(!-3-)/2 RD era 

n ~ < t~ 2w RH/MD eras (1) 

[ exp[— 3H(1 + w)t] inflationary era . 

Recall that w = for matter, while w = — 1 for vac- 
uum energy (cosmological constant) and w = ±1/3 for 
radiation and curvature respectively. 

For concreteness, we shall assume that the individ- 
ual components of the eventual dark matter in our sce- 
nario are described by scalar fields fa with correspond- 
ing masses and widths I\ describing their decays 
into Standard-Model states. For simplicity we shall also 
assume that these widths T, t correspond to processes 
in which a given dark-matter component fa decays di- 
rectly into SM states (i.e., fa — > SM) without pass- 
ing through any other dark-matter components as in- 
termediate states. In other words, we shall assume that 
extra-ensemble decays fa ~~ ^ SM (with combined total 
width Ti) dominate over all possible mira-ensemble de- 
cays fa — > fa + .... It turns out that this is a valid as- 
sumption for many realistic scenarios to be discussed in 
this paper and in Ref. @, and in the Appendix we shall 
discuss what happens when this assumption is relaxed. 

In a Friedmann-Robertson- Walker (FRW) cosmology 
parametrized by a Hubble parameter H(t) ~ t~ x in 
which we assume that our dark-matter component fields 
fa have negligible spatial variations as well as negligible 
self-interactions, these fields will time-evolve according 
to the differential equation 

fa + [3H(t)+Ti]fa + m 2 fa = 0. (2) 

This is the equation for a damped harmonic oscillator, 
with critical damping occurring for 3H(t) + Ti = 2mi, 
note that the single-derivative "friction" term in this 
equation receives two separate contributions, one arising 
from the cosmological Hubble expansion and the other 
arising from the intrinsic decay of fa. As a result, at early 
times for which 3H(t) + Ti > 2m,i, the field fa is over- 
damped: it does not oscillate, and consequently its en- 
ergy density behaves like vacuum energy (with w = — 1). 
By contrast, at later times for which 3H(t) + Ti < 2m,i, 
the field is underdamped: it therefore oscillates, and con- 
sequently its energy density scales appropriately for mat- 
ter (with w = 0). The condition 3H(t) + Ti = 2mi thus 
describes the "turn-on" time at which oscillatory behav- 
ior begins and the field begins to act as true matter. 

Given these observations, we can now sketch how each 
of the abundances for each component fa will behave 



The results in Eq. JTJ follow from the general facts that p ~ 
R~ 3 ( 1 +™) and p cr it ~ H 2 , where R and H are respectively 
the scale factor and Hubble parameter. Recall that these lat- 
ter quantities have the scaling behaviors (R, H) ~ (i 1//2 ,t — 1 ) 
in an RD era; (R, H) ~ (t 2 / 3 ,t _1 ) in RH or MD eras; and 
(R, H) ~ (e Ht , constant) in an inflationary era. 



in our scenario. For concreteness, we shall assume that 
these abundances are all initially established at a com- 
mon time to- Moreover, we shall assume that each com- 
ponent has an initial abundance which decreases as a 
function of its mass. While not all production mecha- 
nisms have this property, we shall see in Sect. Ill that 
misalignment production in particular can accomplish 
this task. 

Immediately upon establishment of these abundances, 
the states in our ensemble can be separated into two 
groups. Those heavier states with masses 3H(to) + 
Ti < 2rrii will all begin oscillating simultaneously. In 
other words, they experience a simultaneous, instan- 
taneous turn-on. By contrast, the lighter states with 
3H(to) + Ti > 2m,i will experience a step- by-step "stag- 
gered" turn-on, with lighter and lighter states crossing 
the turn-on threshold at later and later times. Indeed, 
if we approximate H(t) ss n/3t where k is a constant 
within each era, then a given mode with mass ro, will 
turn on at a time ti w «/2mj. Thus the lightest states 
are necessarily the last to turn on. Indeed, n = 2 for the 
RH and MD eras, while k = 3/2 for the RD era. 

Finally, once these states are all "turned on" and be- 
have as matter with w = 0, their abundances fli will 
evolve as discussed above until such times t ~ = F^ 1 
as these states decay. Specifically, the abundance f2j ~ 
Pi/H 2 associated with each component fa will evolve ac- 
cording to the differential equation 

Oi+ ^3ff + 2^+riJ = . (3) 

[This equation follows directly from Eq. ([2]) upon use of 
the general expression pi — ^(mffa 2 + fa 2 ) and the virial 
theorem ruffa 2 = fa 2 , the latter holding for oscillations 
whose frequencies are large compared with Note 
that if we ignore the decays of these particles (i.e., if 
we set Ti — > 0), the solutions to this differential equation 
are nothing but the results given in Eq. (fTJ) for w = 0. 
Upon decay, however, the corresponding abundance fli 
drops rapidly to zero; this occurs when t ~ Ti, and in- 
deed Ti > 3H + 2H/H at t ~ 75. For simplicity, in the 
rest of this paper we shall approximate such decays as 
occurring promptly and completely at t = r^. However 
this approximation will not be critical for any of our con- 
clusions, and can easily be discarded if needed. 

Combining all these features, we can then sketch the 
salient features of our scenario as in Fig. [TJ In this plot, 
we have taken the time at which the initial abundances 
are established to be during the inflationary era, but 
other time intervals are also possible. We have also as- 
sumed that r.; -c 3H(ti) when 3H(ti) = 2nii, so that the 
decay widths Ti affect the final decay times Tj but not 
the staggered turn-on times U . 

Several features which are clear from this sketch help to 
define and characterize this scenario. First, we see that at 
the present time, there continue to be a plethora of dark- 
matter states. Although each of these has an extremely 
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FIG. 1: A sketch of our dynamical dark- matter scenario in which the dark matter of the universe comprises a vast ensemble of 
individual components with different masses, abundances, and lifetimes. This plot illustrates the evolution of the abundance 
of each dark-matter component as a function of time, assuming that all abundances are initially established at a common 
time (chosen here to be prior to or during the inflationary era), with values that decrease as a function of the component 
mass. For all subsequent times, these abundances scale as vacuum energy until 3H(t) — 2rru, after which point they scale as 
matter. Open circles indicate states which inflate away, while closed circles indicate states which decay into SM particles with 
associated lifetimes that decrease with increasing mass. In our scenario, the lifetimes and abundances are balanced against each 
other in such a way that there continue to exist a plethora of dark-matter states which survive at the present time: although 
each such state has an extremely small abundance Q,i <C 1, they collectively reproduce JIcdm- Nevertheless, because of their 
extremely small abundances, states which have already decayed into SM particles leave negligible imprints on the CMB and 
other observable astrophysical and cosmological backgrounds. An important feature of this scenario is that it is fully dynamical, 
with the composition and properties of the dark matter continuing to experience a non-trivial time evolution before, during, 
and even after the current epoch. 



small abundance (exponentially suppressed on the plot 
in Fig. Q]), they can collectively produce a sizable, 0(1) 
abundance which we choose to identify with the observed 
dark-matter abundance Ocdm ~ 0.23 Q. 

Second, we observe that our dark-matter states are not 
governed by the notion of stability: while some are indeed 
more stable than others, decays of dark-matter states can 
occur throughout the evolution of the universe. This is 
not in conflict with observational constraints because of 
the extremely small abundances of those states which de- 
cay at critical epochs during the evolution of the universe. 
In other words, as discussed above, lifetimes and abun- 



dances are balanced against each other in this scenario. 

Third, as a consequence of these features, we see that 
nothing is particularly special about the present time in 
this framework. Dark-matter states need not be held 
stable until any special moment, and the current age of 
the universe plays no special role in this scenario. Indeed, 
as evident from the sketch in Fig. [TJ dark-matter states 
decay prior to, during, and after the present era. What 
results, then, is a scenario in which the dark matter is 
not frozen in time at the present era, but continues to 
act as a highly dynamical component of an ever-evolving 
universe. 
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C. Characterizing a given dark-matter 
configuration: fitot, 77, and ui e fi 

In general, there are three quantities which we can use 
in order to characterize the configuration of our dark- 
matter ensemble at any instant in time, and to track its 
subsequent time-evolution. We shall begin by introduc- 
ing these three quantities and discussing the relations 
between them. Then, we shall discuss several qualitative 
aspects of their overall time-evolutions. 

1. Fundamental definitions and relations 

The first quantity we shall define in order to character- 
ize the configuration of our dark-matter ensemble is the 
total dark-matter abundance: 

o tot = n * ■ ( 4 ) 

i 

Note that we should include in this sum the contributions 
from only those components in our ensemble which have 
already "turned on" (i.e., have begun oscillating) and 
which are therefore already behaving as true matter. Re- 
stricting such contributions in this way ensures that J7 to t 
can truly be associated with a total dark-matter abun- 
dance. Of course, all of the components in our ensemble 
will eventually "turn on" and behave as dark matter af- 
ter enough time has passed. It is for this reason that we 
shall continue to refer to each component of our ensemble 
as a dark-matter component, regardless of its particular 
turn-on time. 

The quantity £l to t describes the total dark-matter 
abundance at any instant of time. However, we may also 
define a complementary quantity 77 which describes how 
this total abundance is distributed across the different 
components: 



Here fio = maxi{f2i} is defined to be the largest individ- 
ual dark-matter abundance from amongst our ensemble 
of dark-matter states. Thus, 77 measures what fraction 
of the total abundance fitot is not carried by a single 
dominant component. We sec from its definition that 77 
varies within the range < 77 < 1: values of 77 near zero 
signify the traditional situation in which the total dark- 
matter abundance is predominantly carried by a single 
state, while larger values of 77 signify departures from 
this traditional configuration. In this sense, then, 77 can 
also be viewed as quantifying the degree to which our 
scenario deviates from the more traditional dark-matter 
framework at any instant in time. 

Recall that we have assumed for this discussion that 
the more massive components of our ensemble have 
smaller initial abundances, and vice versa. Indeed, this 
assumption is already reflected in the sketch in Fig. [1] 



It then follows that the largest abundance fio in Eq. (|5j) 
will correspond to the lightest component in our ensem- 
ble. However, in the event of a "staggered" turn-on, the 
lightest components in the ensemble are necessarily the 
last to turn on. Indeed, prior to their turn-on times, the 
abundances of these lightest states contribute to the to- 
tal dark- energy abundance rather than to the total dark- 
matter abundance. We must therefore be careful to iden- 
tify Slo as the abundance associated with the lightest of 
those components which have already turned on. As a 
result, even the identity of the component whose abun- 
dance is to be identified with VLq can occasionally be time- 
dependent. 

While both f2 tot and 77 characterize the configuration of 
our dark-matter ensemble at a given instant in time, one 
of the critical features of our dynamical dark-matter sce- 
nario is precisely that it is dynamical — i.e., that these 
quantities have non-trivial time-evolutions. Of course, 
some of this time-evolution is common to all dark-matter 
scenarios, arising due to the Hubble expansion of the 
universe during its reheating, radiation-dominated, or 
matter-dominated eras. There are, however, additional 
time-dependent effects which are unique to our dynam- 
ical dark-matter scenario. For example, one such effect 
dominates the physics of the final, matter-dominated era, 
and arises because the total dark-matter abundance in 
our scenario has been distributed across many individual 
dark-matter components, each with a potentially differ- 
ent lifetime. This phenomenon leads to a slowly falling 
fitot a t late times. Clearly the time-evolution of Q-tot 
during this period is extremely sensitive to not only the 
distribution of the abundances Qi across the different 
dark-matter components within the ensemble, but also 
the decay widths F^ which govern the times at which 
these different components decay. 

For this reason, it will be useful to define a quantity 
which can meaningfully characterize the aggregate time- 
evolution of our ensemble. Moreover, we would like this 
quantity to characterize this time-evolution regardless of 
the particular cosmological era under study, and likewise 
to quantify the extent to which this time-evolution in- 
trinsically differs from that normally associated with the 
cosmological expansion of the universe. 

We have already seen in Eq. ([T|) that the time- 
dependence of the abundance f2i associated with a sin- 
gle dark-matter component can be parametrized by its 
equation-of-state parameter w. It is therefore natural 
to ask what "effective" equation-of-state parameter w e g 
might collectively describe our entire dynamical ensemble 
of dark-matter components. For example, even though 
each individual dark-matter component behaves as mat- 
ter (with w = 0), the decays of these components at late 
times (or even their staggered turn-ons at early times) 
might conspire to produce an effective w-vslue for the 
entire ensemble which is non-zero. Such behavior is il- 
lustrated in Fig. [5] for decays that occur during the final 
matter-dominated era. In all cases and in all cosmolog- 
ical eras, the presence of an effective 7/j e ff which differs 
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FIG. 2: A sketch of the total dark-matter abundance in our scenario during the final, matter-dominated era. Even though 
each dark-matter component individually has w = 0, the spectrum of lifetimes and abundances of these components conspire 
to produce a time-dependent total dark-matter abundance fitot which corresponds to an effective equation of state with w > 0. 



(however slightly) from zero would then signify a depar- 
ture from the traditional dark-matter scenarios. 

We can also understand this at a mathematical level. 
The fact that each individual dark-matter component has 
an abundance which follows the behavior in Eq. (|TJ) with 
w = does not guarantee that their sum fitot must follow 
the same behavior. Indeed, the two effects which can al- 
ter the time-evolution of the sum fitot in our scenario are 
a possible staggered turn-on at early times, and the indi- 
vidually decaying dark-matter components at late times. 
Thus the time-dependence of fitot need not necessarily 
follow Eq. (JTJ with w = 0. 

One possibility, of course, is that fitot will continue 
to follow Eq. ([T]). but with some other effective value 
w e f[. However, even this outcome requires that our indi- 
vidual dark-matter components exhibit certain relation- 
ships between their abundances and lifetimes which need 
not actually hold for our dark-matter ensemble. There- 
fore, in general, we expect that fitot might exhibit a 
time-dependence which does not resemble that given in 
Eq. |TJ) for any constant w e g. Or, to phrase this some- 
what differently, we expect that in general, our effective 
equation-of-state parameter w e ff might itself be time- 
dependent. We therefore seek to define a, function w c g(t) 
which parametrizes a time-dependent equation of state 
for our dynamical dark-matter ensemble as a whole. 

In order to define such an effective function w Q g(t), let 



us first recall that the traditional parameter w is funda- 
mentally defined through the relation p — wp where p 
and p are respectively the pressure and energy density of 
the "fluid" in question. Of course, in an FRW universe 
with radius R, the conservation law for energy density 
dE = —pdV yields the relation d(R 3 p) = —pd(R 3 ), from 
which it immediately follows that (p + p)dR 3 /R 3 = —dp 
or 3(p + p)d\ogR — —dp. Recognizing p + p = (1 + w)p 
and dlogR — Hdt where H is the Hubble parameter, we 
thus have 



3H(l + w) = - 



d log p 
dt 



(6) 



This is a completely general relation which makes no as- 
sumptions about the time- (in) dependence of w. We may 
therefore take this to be our fundamental definition for 
w e s(t) — i.e., 



W e s(t) 



( 1 dlogptot 
\3H dt 



d log fitot 
d log t 

d log fitot 
d log t 



1 



for RH/MD eras 



for RD era 



(7) 



Note that while our derivation has thus far been com- 
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pletely general, we have specialized to specific cosmolog- 
ical eras in passing to the final expressions in Eq. (J7J. 
Specifically, we have written p tot — O to tPcrit and taken 
3H ~ n/t where n = 2 (RH/MD), « = 3/2 (RD). 

The final expressions in Eq. ([JJ are easy to interpret 
physically, since the double-logarithmic derivatives which 
appear in these expressions are nothing but the slopes in 
the sketches in Figs. []] and [5] However, the important 
point of this derivation has been to demonstrate that w e ff 
defined as in Eq. continues to have a direct interpreta- 
tion as a true equation-of-state parameter relating energy 
density and pressure, even when w c s is time-dependent. 
No other definition of w e s would have had this property. 

The results in Eq. ([7]) provide a relation between w e s 
and fitot(i)' However, it is also possible to derive a sim- 
ilar relation between w e ff and 77. Assuming that we re- 
strict our attention to periods of time after all staggered 
turn-ons are complete (so that the identity of the dark- 
matter component associated with Oo is fixed), it triv- 
ially follows from the definition of 77 in Eq. © that 



dlog(l - 7/) 
dlogt 



d log Otot 

dlogt 
d log O to t 



RH/MD eras 



RD era 



dlogt J 2 
Using the results in Eq. ([7]), we therefore find that 



w e s(t) 



dlog(l - 77) 



dlogt 
dlog(l - 77) 



dlogt 



RH/MD eras 
RD era . 



(8) 



(9) 



It therefore follows that decreasing 77 corresponds to pos- 
itive Weff, and vice versa. As a self-consistency check, we 
also observe that the standard paradigm — which has 
7] = for all times — has w e g = for all times as well. 



2. fltot, r\, andw e s: Qualitative time- dependent behaviors 

Having defined the three quantities {Otot , V> w ob} that 
we will use in order to characterize a general configura- 
tion of our dark-matter ensemble, we now seek to under- 
stand the time-dependence that will be generically exhib- 
ited by these quantities across the different cosmological 
eras. As we shall see, several qualitative observations 
can be made even without further assumptions concern- 
ing the individual abundances and decay widths within 
our ensemble. 

Given the sketch in Fig. Q] it is perhaps easiest to un- 
derstand the qualitative behavior of O to t as a function 
of time. There is, of course, the time-dependence for 
Otot which can be associated with the regular Hubble 
expansion of the universe and which causes fi to t to in- 
crease during the radiation-dominated era. This aspect 
of the time-dependence is generally common to all dark- 
matter scenarios. However, as discussed above, there 



are two additional effects which are specifically associ- 
ated with our dynamical dark-matter scenario and which 
cause fitot to experience a further time-dependence. The 
first of these is the possibility of a "staggered" turn-on 
across the different dark-matter components in our en- 
semble — i.e., the possibility that some components will 
remain longer than others in a state in which their abun- 
dances Hi grow rapidly as ~ t 2 but contribute to the total 
dark-energy abundance rather than to Otot- This feature, 
when present, will thus generally cause fi tot to experience 
a more gradual growth (but also a greater eventual max- 
imum value) than would otherwise occur in a more tradi- 
tional dark-matter scenario. However, at later stages of 
the cosmological evolution (particularly during the final, 
matter-dominated era), we see that our different dark- 
matter components have a broad spectrum of lifetimes 
and decay widths. This causes Otot to experience a slow 
step- wise decline before finally reaching zero upon the de- 
cay of the last-surviving dark-matter component within 
the ensemble. 

Similar qualitative arguments also apply to the time- 
evolution of 77. An initial value of 77 is implicitly deter- 
mined once the abundances for each of the dark-matter 
components are established. Of course, if this occurs dur- 
ing an inflationary period, it is possible that certain more 
massive components will have inflated away by the time 
the inflationary period ends. If this is the case, then we 
may regard the "initial" value of 77 to be the value of 77 
at the end of the inflationary period. 

In general, after that point, the evolution of 77 can ex- 
perience as many as three distinct phases. Let t\ denote 
the time at which the last (lightest) dark-matter compo- 
nent has turned on, and let ti denote the time at which 
the most massive dark-matter component decays. As- 
suming t<i > ti, there are therefore three distinct time 
intervals which become relevant. 

During ti < t < t%, each of the individual dark- 
matter abundances experiences a common overall time- 
dependent scaling behavior as the universe evolves. As 
a result, the ratio between the abundances of these com- 
ponents remains fixed. In other words, 77 remains frozen 
during this period (even though O to t may continue to 
vary). 

For t > t%, by contrast, it is clear that the decays of the 
more massive dark-matter components have the cumula- 
tive effect of decreasing O to t without altering Hq. What 
results, then, is a step-by-step, threshold-by-threshold 
decline in the value of 77. This process continues un- 
til only a single dark-matter component survives and 77 
reaches zero. By assumption, in our scenario this will not 
happen until the distant future, at t = T^ 1 t now where 
Ti is the decay width of the second-lightest dark-matter 
component. 

However, during the period t < t\, a "staggered" turn- 
on for the individual dark-matter components can also 
generally induce a non-trivial time-evolution for 77. In- 
deed, each time a new lightest dark-matter component 
turns on, its abundance Og suddenly contributes to the 
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total dark- matter abundance f2tot- This abundance fl' 
also displaces the previous largest individual abundance 
f?o- We therefore find that with each such successive 
turn-on, rj experiences a shift in its value: 



77=1- 



rf = 1 - 



•nf, 



(10) 



By assumption, f2g > f2o- However, it is easy to see 
from Eq. (ITU)) that rj' < r\ only if il' > Qo/rj. Since 
< ?7 < 1, we see that this condition is guaranteed to be 
satisfied only if rj = 1, and guaranteed not to be satisfied 
only if rj = 0. In all other cases, this condition may or 
may not be satisfied, and this will cause rj to either de- 
crease or increase, respectively. We also observe that in 
a very rough sense, rj tends to stabilize and avoid cither 
the rj = 1 or the rj = extremes: as rj — > 1, it becomes 
easier and easier to satisfy the constraint that drives rj 
lower, while as rj — ¥ 0, it becomes easier and easier to 
satisfy the constraint that pushes rj higher. Indeed, if 
we imagine that each new abundance Q' somehow has a 
random value greater than the previous Oo, we can envi- 
sion an "oscillatory" behavior in which rj varies between 
its two limits. Unfortunately, we cannot be more specific 
about this behavior without knowing something further 
about the individual abundances that exist during such 
a staggered turn-on phase. 

We therefore conclude that rj will take an initial value 
once the abundances are established, and that this value 
can then undergo a non-trivial time dependence if there is 
an initial period during which a staggered turn-on occurs. 
After the staggered turn-on is complete, rj will remain 
frozen until late times when our individual dark-matter 
components begin to decay. This will then cause rj to fall 
monotonically, ultimately vanishing when only the single 
longest-lived dark-matter component remains. However, 
it is regarded to be a fundamental property of our sce- 
nario that rj is nevertheless presumed to be significantly 
different from zero at the present time. 

Finally, we turn to the behavior of w e g as a function 
of time. However, given the relations in Eq. @, it is 
relatively straightforward to map out the rough time- 
dependence of w e g. During the staggered turn-on phase 
(i.e., t < t\), we have seen that rj may either increase 
or decrease; this implies that w e s may be either negative 
or positive. Moreover, the fact that rj tends to stabilize 
during this period, avoiding its extreme values at rj = 
or rj = 1, implies that w e g will likewise tend to stabilize 
around zero, with positive values of u; c ff ultimately fol- 
lowed by negative values, and vice versa. As indicated 
above, however, this assumes that each new abundance 
51q that turns on has a value which is greater than the 
previous fio but is otherwise somewhat random. 

During the period t\ < t < i 2 , by contrast, the be- 
havior of ui e ff is far simpler to describe: we simply have 
w c fF = 0. This is completely in accord with our obser- 
vation that rj stays frozen during this period, and that 
our dynamical ensemble is behaving as ordinary dark- 
matter during this period (except with a non-zero value 



of rj). Finally, for the period t > t 2 after our individual 
dark-matter components have begun to decay, we have 
argued that rj is monotonically decreasing. This implies 
that w e e is strictly positive during this period, a feature 
which is illustrated in Fig. [2] and which again serves as 
a cosmological "smoking gun" for our dynamical dark- 
matter scenario. 



D. A signature of dynamical dark matter: 
Time-evolution of Qtot, rj, and w a g during the final 
matter-dominated epoch 

As just discussed, one of the most important signatures 
of our dynamical dark-matter framework is the fact that 
the total dark-matter abundance 51 to t is a time-evolving 
quantity — even during the current matter- dominated 
epoch. Within such a framework, it is therefore only to 
be regarded as an accident that this quantity happens to 
reproduce a specific observed value ficDM ~ 0.23 at the 
present time. 

With only a few additional assumptions, it turns out 
that we can explicitly calculate the time-evolution of the 
total dark-matter abundance fitot during this epoch. We 
can also explicitly calculate the time-dependence of 77, 
and the resulting equation of state w c s(t). In the rest 
of this section, we shall therefore concentrate on the fi- 
nal matter-dominated epoch. Indeed, this is the epoch 
during which a non-trivial time-evolution for f2 t ot arises 
only because of the decays of the individual dark-matter 
components within our ensemble. 

Within this era, each dark matter component 4>i has 
a relative abundance fli which remains constant until it 
decays at a time t ~ Tj = iy 1 . Taking this decay to be 
nearly instantaneous, we can thus write 



(11) 



whereupon we see that 
dtttot (t) s-^ „ d 



dt 



-t) 



■^QiSfc-t) (12) 



where we have defined il to t(t) = an( ^ use( i the 

relation dO(x)/dx = 8(x) where S(x) is the Dirac 5- 
function. In the limit that we truly have a large number 
of dark-matter states, we can imagine that the spectra 
of decay widths T.; and decay times Ti = 1 are nearly 
continuous, with continuous variables T and r. With this 
approximation, we can view fij as a continuous function 
£1(t) and convert the sum over states to an integral, i.e., 



E 



dr n T (r) 



(13) 



where n T (r) is the density of dark-matter states per unit 
of t, expressed as a function of r. Eq. (TT2"j) then becomes 



dCl tot (t) 
dt 



(frft(-r) n T {T)5(r-t) 



(14) 
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In general, the quantities n[r) and f2(r) are unspeci- 
fied, their properties depending on the particular dark- 
matter scenario under study and the specific features of 
our dark-matter ensemble. However, it will prove con- 
venient to parametrize these quantities in terms of their 
scaling behaviors as functions of T: 



fi(r) w AY a , 



n r (r) 



BY 



/3 



(15) 



with overall (generally dimensionful) coefficients (A, B) 
and scaling exponents (a,/?). Since the abundances of 
states in our scenario generally have an inverse relation 
to their decay widths, we expect that a < 0. Note that 
np in Eq. f| 15[) is the density of states per unit of T, 
whereupon it follows that 



"r 



dT 



= T 2 n v 



(16) 



We thus find that fi(r)n T (r) - ABT a+l3+2 , or equiva- 
lently 0(r)n r (r) - ABr~ a ~P~ 2 . Use of Eq. (0 then 
leads to the result 



dtttotjt) 
dt 



ABt 



-a— 0-2 



(17) 



Imposing the condition that f2 to t = ^cdm at the present 
time t = t now and assuming that a + /3 =/= —1 then leads 
to the solution 



^tot(i) — ^cdm + 



AB 



a + P + l 



-a—fi-l _ f-u-p-1 



For a + j3 = — 1, by contrast, we have the solution 



fitot(t) = ^CDM-ABlog 



tn 



(18) 



(19) 



Under the assumptions in Eq. (jTSJ), the results in 
Eqs. (fTE)) and (IT9")) are completely general in a matter- 
dominated era. Moreover, it is clear from Eqs. (fT8f and 
(fP9|) that in all cases, fitot decreases with time. This is 
precisely as expected, since all of the time dependence of 
f2tot in a matter-dominated era arises due to the decays 
of the individual dark-matter components within the en- 
semble. Notice that some of these functional forms for 
fitot actually predict that Q to t(t) will eventually hit zero; 
this is also not unexpected, since this corresponds to the 
final decay of the last remaining dark-matter component 
in the ensemble. Needless to say, we should not consider 
any of these f2tot(£) functions beyond the times when they 
might hit zero. Nevertheless, as long as our dark-matter 
ensemble obeys the scaling laws in Eq. (fT"5|) , the functions 
given in Eqs. (fT8"|) and (fT9"|) correctly describe the behav- 
ior of the corresponding total dark-matter abundances 

ntot(*). 

Given the results in Eqs. (IT8l) and (fl~9| as well as the 
definition in Eq. ([7]) for a matter-dominated era, we can 
also obtain a solution for the time-dependent equation- 
of-state parameter w c g(t) associated with our ensemble 



of decaying dark- matter states. For x = a + /3 ^ — 1, we 
find 



w e s(t) 



(1 + x)w* 



2w* + (1 + x - 2iu*)(i/inow) 1+x 



(20) 



where 



AB . . 

w* = w cB (t aow ) = — -j^ . (21) 

^"CDMtnow 

Note that for w» 1, this result is fairly well- 

approximated by 



t 



W e s{t) ~ W* 



By contrast, for x — — 1, we instead obtain 



1 - 2w* log(t/t 

now I 



where 



w* = w cS (t now ) 



AB 



2VL 



CDM 



(22) 



(23) 



(24) 



The behavior of the results in Eqs. and (HH de- 
pends critically on the relationship between x and u>*. 
For 1 + x < 2w*, we find that u> c ff always increases mono- 
tonically as a function oft before reaching w„ at t — t aow . 
By contrast, for 1 + x > 2w*, this function decreases 
monotonically before reaching tu* at t = t now . Finally, for 
1+x = 2w„, we have the exact result that w e s(t) = u>* for 
all t. This (admittedly fine-tuned) case illustrates that 
it is possible to achieve a time-independent equation- 
of-state parameter w e f{ = under the assumptions in 
Eq. (|15l) , and moreover that this value of can be tuned 
to any positive value desired. This is indeed the situa- 
tion illustrated in Fig. [21 which is plotted for a < and 
P > 0. 

The above qualitative descriptions indicate the history 
of w e ff(i) prior to the present day. However, in general, 
this same increasing or decreasing behavior continues for 
t > tnow (i.e., through and beyond the current epoch), 
with one important caveat: for 1 + x < 2w*, we see 
that w c ff(t) not only continues to increase, but eventu- 
ally hits a pole. However, such poles represent the lo- 
cations at which the corresponding f2 to t-functions have 
zeroes. These poles are therefore unphysical, signalling 
the decay of the last component within our dynamical 
dark-matter ensemble, and we can restrict our analysis 
of these functions to times preceding these critical values. 

If our dynamical dark-matter scenario is to be in 
rough agreement with cosmological observations, we ex- 
pect that ui* today should be fairly small (since tradi- 
tional dark "matter" has w = 0). We also expect that 
the function w e g(t) should not have experienced strong 
variations within the recent past. This suggests that sit- 
uations with x < — 1 are likely to be phenomenologically 
preferred over those with x > — 1, since having x < — 1 
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ensures that < w Q g(t) < tw* for all t < t now . Indeed, the 
more negative x becomes, the closer to vanishing w e g(t) 
remains before finally reaching w* at t = t now . However, 
depending on the detailed properties of the particular 
realization of our dynamical dark-matter scenario under 
study, values of x near —1 or slightly above may also be 
phenomenologically acceptable. 

Finally, we may also use these results to solve for 77 as 
a function of time. For x = a -\- /3 ^ — 1, we find 



2w* + [ry*(l + x) - 2w*](t/t now ) 



l+x 



2 Wsr + [1 + x - 2w*](t/t now ) 1 +* 
where is given in Eq. (1211) and where 



v* = v(t 

Likewise, for x ■ 

v(t) = 



= 1 



fin 



ficDM 

1, we have 

r]* - 2w* log(f/f n ow) 

1 - 2w* log(f/f now ) 



(25) 



(26) 



(27) 



where w* is given in Eq. (J24J). It is not surprising that 
r), unlike w e g, depends on two independent dimension- 
less quantities u>* and 77*, since the very definition for 77 
introduces a new quantity fin which had not previously 
appeared. 

Note that all time-dependence for 77(f) cancels, with 
rj(t) m 77* for all x, if either w* — > or 77* — > 1. This 
makes sense, since in the first case fitot does not change 
while in the second case fin ~> 0- These are the only 
two ways in which 77 can remain constant. In all other 
cases, however, 77(f) is always a decreasing function of 
time, as expected. We also see from Eqs. ([25]) and (|27|) 
that 77(f) — > 1 for all x > — 1 as t/t now — > 0. Indeed, this 
holds regardless of the values of w* or 77* . 

While these characteristics successfully conform to our 
expectations concerning the behavior of 77(f), there are 
some features that the functional forms in Eqs. (|25[) and 
(|27p do not accurately capture. For example, if x + 1 < 
2w*/?7*, these functions predict that 77(f) will eventually 
become negative beyond a certain late time. Moreover, 
while these functions resemble those for w e g in that they 
properly capture the pole that results when fi to t — > 0, 
they do not necessarily approach 77 — > before hitting 
this pole. 

The reason behind these failures is easy to understand. 
Unlike fitot and w e g, the quantity 77 has a special charac- 
teristic not shared by the others: it is sensitive not only 
to fitot, but also to fin. Features such as having 77 reach 
zero but not become negative are extremely sensitive to 
the value of fin and the fact that fitot must exactly hit 
fio after all but the lightest dark-matter component have 
decayed. Indeed, these features are extremely sensitive to 
the fine-tuned and ultimately discrete nature of the light- 
est dark-matter components, and this is precisely the sort 
of information that our scaling assumptions in Eq. (|15[) 
are incapable of modelling. Thus, while we may view the 



functions in Eqs. (|25|) and (|27|) as being reasonably ac- 
curate for most portions of the cosmological evolution in 
the matter-dominated era, we should not maintain this 
expectation beyond a certain time when all but a few 
dark-matter components have decayed. 

To summarize, then, in this section we have presented 
a dynamical multi-component dark-matter scenario in 
which individual component abundances and lifetimes 
are balanced and distributed across the components in 
such a way that constraints from BBN and other back- 
grounds are potentially satisfied. An important part of 
this scenario is the proposition that both 77 and w e g are 
different from zero at the present time, the former signif- 
icantly so, and that components of the dark matter are 
actively decaying prior to, during, and beyond the cur- 
rent epoch. As a result, cosmological quantities such as 
fitot experience a time-evolution which transcends that 
due to the ordinary expansion of the universe. 



III. DYNAMICAL DARK MATTER MEETS 
THE INCREDIBLE BULK 

Thus far, we have done little more than present a new 
framework for dark-matter physics. In particular, we 
have not yet demonstrated that an ensemble of dark- 
matter states can easily be assembled in which the in- 
dividual component abundances are naturally balanced 
against lifetimes in a well-motivated way. In this section, 
however, we shall demonstrate that an infinite tower of 
Kaluza-Klein states propagating in the bulk of large extra 
spacetime dimensions naturally constitutes an ensemble 
of states with the desired properties. As we shall see, 
this occurs because KK towers obey a special "balanc- 
ing" constraint which relates the lifetimes of individual 
KK modes to their abundances. Specifically, we shall 
demonstrate that the KK modes within a generic KK 
tower exhibit abundances fij and SM decay widths Ti 
which obey an inverse relation of the form anticipated in 
Eq. $T5\), i.e., 



fiiH 



constant 



(28) 



for some a < 0. This constraint ultimately emerges as a 
consequence of the non-trivial interplay between physics 
in the bulk and physics on the brane. 



A. General setup 

For simplicity, we shall consider our spacetime to take 
the form M.4, x S1/Z2, where A^4 denotes ordinary four- 
dimensional Minkowski spacetime and S1/Z2 denotes a 
line segment which is realized as a 7Li orbifold of a cir- 
cle of radius R. We shall take z M = (x^,y) to denote 
the coordinates on this spacetime, with the Z2 orbifold 
action identified as y — > —y, and imagine that the Stan- 
dard Model is restricted to a brane at the fixed point 
y = 0. We are therefore considering a "toy" ADD-like 
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scenario [10( with a single flat extra dimension. Despite 
the simplicity of this toy model, we are making no as- 
sumptions at this stage about relevant mass scales or the 
full number of extra spacetime dimensions that might ac- 
tually exist in a more fully realized scenario. Indeed, we 
believe that most of the desired properties that emerge 
from this scenario are likely to be retained if we imagine 
that our spacetime contains additional extra dimensions, 
or is warped @ rather than flat. 

In such a scenario, all fields that propagate in 
the "bulk" are necessarily singlets with respect to all 
Standard-Model gauge forces. As a result, such fields can 
have at most highly suppressed (e.g., gravitational) inter- 
actions with the Standard-Model fields, and thus appear 
as dark-matter candidates. Such fields might include the 
graviton, axion, and other moduli fields. For simplicity, 
we shall consider the case in which the bulk field is a 
five-dimensional scalar but we shall make no further 
assumptions about its properties. 

Neglecting gravity, and with tpi generically denoting 
the Standard-Model fields, we see that such a scenario 
therefore has an action of the form 



S = Jd 4 xdy [£ b uik($ 



) + %)£ branc (^,$)] . (29) 



In general, we may assume that our five-dimensional bulk 
action takes the form 



£buik = \d M $*d 



M 



±M 2 |$| 2 



(30) 



where <9m denotes a five-dimensional derivative and 
where M is an unspecified bulk mass. In certain cases, 
specific symmetries may restrict us to the case with 
AI = 0, but we shall leave M general until further notice. 

Likewise, the brane action will generically consist of 
two contributions — the usual Standard-Model action 
Csm , and an action £j nt which arises due to the interac- 
tions between $ and the Standard-Model fields: 



£brane — -^SM + Mr 



(31) 



In general, there are two types of interactions which will 
concern us. The first class of interactions result in explicit 
couplings between $ and the Standard-Model fields, and 
will ultimately be responsible for allowing $ to decay into 
Standard-Model states. We shall discuss such interac- 
tions in Sect. III.B. There is, however, another possible 
type of interaction term which can also appear within 
£; n t: this is a possible "brane mass" for $ itself, i.e., 



D -±m 2 |$| 



(32) 



Such a brane mass can emerge as an effective operator 
arising due to perturbative or non-perturbative dynamics 
wholly restricted to the brane. Note that this brane-mass 
term must not be confused with the primordial bulk mass 
that appears in Eq. (|30p ; rather, this term has its origins 
within the physics on the brane itself, and appears as 
part of £brane within Eq. (|2"9")l rather than £buik- 



These minimal assumptions are already sufficient to 
permit us to understand the nature of the resulting 
Kaluza-Klein spectrum for $. Indeed, the followin g re - 
sults are similar to those previously obtained in Ref . |11| . 
As appropriate for compactification on the line segment 
S x jlL-i , we can decompose our five-dimensional field $ in 
terms of an infinite tower of four-dimensional modes <f>k , 



V2ttR fc=Q V « 



the normalization factors 



r 1 for k = 
\V2 for k > 



(33) 



(34) 



are designed to ensure that each mode <pk has a canon- 
ically normalized kinetic term in the resulting four- 
dimensional theory. We then find that 

S = J d 4 xdy [\d M $*d M $- ±M 2 |$| 2 - \8(y)m 2 $> 2 } 




DO 

5 E 

k,£=0 



M 2 



hi 



Mi (35) 



where the Kaluza-Klein (mass) 2 matrix is given by 



Ml 



kt 



R 2 



+ M S M + r k r e m 2 . (36) 



Given these results, we see that this mass matrix would 
have been diagonal were it not for the brane mass term. 
This in turn implies that the KK mass eigenstates <f>\ 
necessarily differ from the KK momentum eigenstates 4>k 
— i.e., there is a non-trivial mixing that is induced as a 
result of the KK mass. This mixing turns out to be crit- 
ical for our analysis. In general, we may characterize the 
degree of mixing in terms of the dimensionless parameter 



y 



l 

mR 



(37) 



For y 1 the mass matrix is essentially diagonal; 
this is what trivially occurs, for example, in the four- 
dimensional R — > limit in which the excited KK modes 
decouple. By contrast, in the opposite limit y <C 1, 
the mixing is essentially maximal across all of the eigen- 
modes. 

It is possible to describe the solutions for the eigenval- 
ues A 2 of the (mass) 2 matrix in closed form, and thereby 
obtain explicit expressions for the corresponding mass 
eigenstates. The eigenvalues turn out to be the solutions 
to the transcendental equation 



vrm 2 i?cot (vri?VA 2 - Af 2 ) = \J A 2 - M 2 



(38) 



If m were zero (i.e., no brane mass), the solutions to this 
equation would be nothing but the expected eigenvalues 



A 



M 



n 

R 2 



(39) 
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and more generally this remains approximately true when 
m <C 1/R, i.e., when y ^> 1. What is perhaps surprising, 
however, is that the presence of a non-zero brane mass 
does not result in a further additive shift in this mass 
spectrum for the KK tower (as does the bulk mass term) , 
but instead distorts the lower mass eigenstates in the 
tower so that they approximately follow the alternate 
spectrum 



A 2 = M 



R 2 



(40) 



Remarkably, this is precisely the spectrum which we 
would normally associate with a five-dimensional field 
$ which is taken to be anti-periodic (rather than peri- 
odic) around the extra-dimensional circle prior to orb- 
ifolding! Indeed, for general values of y, the solutions A n 
of Eq. ([35]) tend to follow the spectrum in Eq. (g0) for 
n <C 7r/y 2 , while they follow the spectrum in Eq. (|39|) 
for n 3> Tr/y 2 and smoothly transition between the two 
spectra for intermediate values n ~ ir/y 2 . As we have dis- 
cussed above, this unusual behavior is the consequence of 
the non-trivial interplay between brane and bulk physics, 
and may have applications beyond its appearance here. 

For each mass eigenvalue A, we can also solve for the 
corresponding mass eigenstate \<j)\) as a linear combina- 
tion of the KK- momentum eigenstates \<f>k). We find the 
exact result [ill 



\<t> 



°° r X 2 



where we have defined the dimensionless eigenvalues 



(41) 



and where 



-4> 



A = VA 2 - M 2 /r 



V2 



(42) 



(43) 



1 + TT 2 /y 2 + X 2 



Given these results, it is straightforward to convert 
between the mass-eigenstate basis \4>\) and the KK- 
momentum basis \4>k)- It turns out that there are two 
specific groups of matrix elements involved in this con- 
version which will be of particular interest to us. The 
first involves the KK zero-mode 4>k=o , for which we have 
the matrix elements 



1.5 




A/y 



FIG. 3: Values of Ax (falling curves) and X 2 A\ (rising curves), 
plotted as functions of the mass eigenvalues X/y = XR for y = 
ty (black), y = ^/n (blue), and y — 1 (red). For each y, there 
are only a discrete set of corresponding allowed eigenvalues 
A (indicated with solid dots); note that the quantity X/y — 
VA 2 — APR takes values closer to 2 + 1/2 near the bottom 
of each tower and shifts to values closer to 2 as A increases. 
In each case, we see that A\ falls with increasing A, while 
X 2 A\ increases with increasing A and ultimately reaches an 



asymptote X 2 A\ — > y/2 as A — > oo. 



For this projection field, we likewise have the matrix el- 
ements 



00 2 \ 2 

k=0 a 



y 



■ cot I — I A x = X 2 A X 



(46) 



Ay 



(44) 



where we have made use of Eq. (|38l) in the final equality. 

In Fig. [3j we plot the values of A\ and X 2 A\ as func- 
tions of A and y. We see that A\ falls with increasing 
A, while X 2 A\ increases with increasing A and ultimately 
reaches an asymptote A 2 ^4a — > V% as A — > oo. Moreover, 
we see that larger and larger values of A are needed to 
reach this asymptote as y decreases. 



However, the second concerns the ■projection of the five- 
dimensional bulk field $(?/) onto the Standard-Model 
brane at y = 0, i.e., 

oo 

^ = <%)Uo = X r ^ ■ ( 45 ) 

fe=0 



B. Balancing V versus f2 

We now address the central feature underpinning dy- 
namical dark matter: the balance between the SM decay 
widths T\ associated with each KK mass eigenstate and 
the corresponding cosmological abundances Q\. As we 
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shall show, an inverse relation of the form anticipated in 
Eq. ([25)1 naturally emerges across the entire KK tower. 



1. Abundances £l\ 

We begin by focusing on the different cosmological 
mode abundances fl\ that can arise in such a scenario. 

During the course of the evolution of the universe, 
there are many production mechanisms through which 
the different KK states might come to be populated and 
thereby acquire non-zero abundances. One such method, 
for example, is thermal production; another relies on 
purely geometric effects (e.g., topological defects such as 
cosmic strings and domain walls) and the decays associ- 
ated with them. However, there is also a third production 
mechanism which exists in cases where the bulk mass M 
happens to vanish: this is so-called "misalignment pro- 
duction" . 

In many string-theoretic contexts, bulk fields often do 
have vanishing bulk masses. Such fields often include 
gravitational and/or geometric moduli fields; they also 
include various axion-like fields. Moreover, as we shall 
demonstrate, the predictions of misalignment production 
are rather straightforward to calculate, and are fairly 
generic for bulk fields as a whole. We shall therefore 
take M = in what follows, and restrict our attention 
to abundances established through misalignment produc- 
tion. 

It is easy to understand the physical underpinnings 
of misalignment production within the framework of dy- 
namical dark matter. Prior to the brane dynamics that 
establishes the brane mass m, the fact that M = implies 
that our theory exhibits a five-dimensional shift symme- 
try $ — > $ + c, where c is a constant. As a result, any 
value for ($) is equally likely to occur: 

(*) = OlT (47) 

where 6 is a random 0(1) dimensionless coefficient and 
where /$ is a mass scale (or decay constant) associated 
with the five-dimensional $ field in the bulk. Decom- 
posed into KK eigenstates via Eq. (|3"3"|) . this non-zero vev 
for the five-dimensional field $ implies a non-zero vev for 
the KK zero mode: 

(<h) = U - (<t>k) = for all k > (48) 

where 

U = y/^RfT ■ (49) 

Note that all of the higher KK modes tfik with k > must 
have vanishing vevs as a result of the five-dimensional 
shift symmetry. 

This is the situation that exists prior to (i.e., at en- 
ergies higher than those associated with) the brane dy- 
namics that establishes the brane mass m. However, once 



this brane mass is established, we must shift to the mass- 
eigenstate basis, whereupon we see from Eq. (1441) that 
Eq. P5|) now takes the form 

(4>x) = OAxU for all A . (50) 

Thus, we see that all of the mass eigenstates will gener- 
ally have non-zero values. Of course, the fact that that 
these vevs are all related through the coefficients A\ is a 
reflection of our original five-dimensional shift symmetry 
in the bulk. 

The dynamics that establishes the brane mass m also 
leads to a non-zero energy density p associated with 
the configuration in Eq. ([SU]) . In general, the four- 
dimensional energy density p\ associated with each mode 
4>\ is given by p\ — \\ 2 (4>\) 2 . Given Eq. (jSTJ)) , we thus 
have 

Px = \e 2 \ 2 A 2 J 2 . (si) 

Of course, at any moment in the evolution of the universe, 
the critical energy density is given by 

Pcrit = 3M 2 H 2 (52) 

where Mp = (87rGjv) -1 / 2 is the reduced Planck scale and 
where H is the Hubble parameter. The initial abundance 
SI a = px/pcrit associated with the 4>\ mode is thus given 
by 

6 I M P H I v ' 

This is, in fact, a completely general result. 

We shall let t denote the time at which this initial 
abundance is established by the brane dynamics. Thus 
Q\(to) — The next question, however, is to de- 

termine the corresponding value of £l\(t now ). In order 
to do this, we see from Fig. [T] that we must make some 
assumptions about whether to is situated during the re- 
heating, radiation-dominated, or matter-dominated eras, 
and whether the <j)\ mode experiences an instantaneous 
turn-on at to or a staggered turn-on at a time t\ > to- 
There are therefore six different cases to consider. 

For simplicity we shall assume that for the modes <f)\ 
which are part of a staggered turn-on, the corresponding 
turn-on time £a occurs at the threshold 3H(t\) = 2A. 
(Of course, if t\ < to, then such modes turn on only 
at to.) We shall also assume that all modes in a given 
tower actually turn on during the same era as to > so that 
our turn-on "cascade" down the tower does not cross a 
boundary between two different eras. Finally, we shall 
assume that within each era, H(t) takes the approxi- 
mate form H(t) = k/M where K — 2 for the reheat- 
ing and matter-dominated eras and k = 3/2 during the 
radiation-dominated era. This implies that £a = k/2X. 
Note that this approximate form for H(t) is generally 
valid at relatively late times within each era, and we 
shall disregard all 0(1) "threshold" effects associated 
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with the boundaries between different eras. Thus, we 
shall implicitly take Q x to be a continuous function of 
t, as sketched in Fig. [TJ and we shall therefore disregard 
all A- independent 0(1) numerical coefficients in those ex- 
pressions for f^A^now) which follow. 

Clearly, if to occurs during the final matter-dominated 
era (i.e., if to > £mre)j then modes which turn on instan- 
taneously (i.e., modes with t x < to) will have abundances 



and 



^A (^now) 



x x n\ 



(o) ( *a\ f £mre\ 

M / V ^RH / 



1/2 



ftp | ( ^MRE \ 

Mp j I £rh / 



1/2 



(60) 



^A(inow) ~ X x n[ Q) ~ ~^A\X X I Aj (mt ) 2 

where X x denotes the expected damping factor due to 
dark-matter decays: 



(55) 



By contrast, for those modes which experience a stag- 
gered turn-on (i.e., modes with t x > to), this result be- 
comes 



^A(inow) ~ X x n[ Q) fe) ~ AlX x l^-\ (56) 



where we have substituted the result t\ ~ 1/A in passing 
to the final expression. 

By contrast, for i within the radiation-dominated era 
(i.e., £rh ^ to < ^mre), these two cases are instead given 
by 



^A^now) ~ X\ 0^ \~Y~J 



\ 2 A\X X ^Aj (^ ) 3 / 2 Mmre) 1/2 

(57) 



Interestingly, of all six cases, this is the only one which 
yields a result for f^A^now) which is parametrically inde- 
pendent of the scale m. 

It is also instructive to examine the manner in which 
these results scale with A. Surveying Eqs. (j5"4")) through 
(f6"0"|) , we see that the dependence of Vl x on A follows only 
three different patterns, depending on the specific turn- 
on behavior experienced by the KK mode in question and 
the era during which it takes place: 



\ 2 A\ 
A 1 / 2 A? 



instantaneous 
staggered (RD era) 
staggered (RH/MD eras) 



(61) 



Under the assumption of misalignment production, this 
result is exact and completely general. However, given 
the definition in Eq. (|43p. we may approximate 



Ay 



f 1/A 
1 1/A 2 



for A < y/l + n 2 /y 2 



for A > y/l + ir 2 /y 2 



(62) 



In the future, we shall refer to these two approximation 
regimes as the small-A and large-A regimes; note that 
while there always exists a large-A regime, the existence 
of a small-A regime depends on the value of y. We then 
find that the results in Eq. (jBTj) lead to the large-A scaling 
behaviors 



A" 



instantaneous 



and 



A 7 / 2 staggered (RD era) 

A -4 staggered (RH/MD eras) 



(63) 



to) \ t x 
\^AiX x (A\ (mW) 1/2 



Finally, for t within the reheating era (i.e., t < £rh)j 
these two cases are instead given by 



n x (t now) ~ Ma 0) (^f) 



1/2 



i2 «^(^) <"*■>' (^) 



as well as the small-A behaviors 




instantaneous 

staggered (RD era) (64) 
staggered (RH/MD eras) . 



(58) Indeed, the only A-dependence which is not included in 
these results is that which appears through the decay 
widths in the JTA-factors in Eqs. ([54]) through (|60|) . How- 
ever, these X x factors only express the physics of the 
eventual dark-matter decay processes; they play no role 
in determining the mode abundances that exist prior to 
decay, which is our main interest in this discussion. As 
a result, we shall disregard these X x factors in what fol- 

1/2 lows, understanding that our analysis is primarily ap- 

propriate for the period that exists prior to the onset 
of decays of the KK states with significant abundances. 

(59) (Indeed, the final period of KK decays will be discussed 
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in Sect. III.B.3 after we have analyzed the behavior of 
the decay widths T\ in Sect. III.B.2.) 

Given the individual abundances il\ in Eqs. (|54|) 
through (j6"0|) , we can now calculate the values of both 
f2tot and 77 that exist prior to the onset of significant KK 
decays. Recall that Qtot is nothing but the sum over all 
of the individual abundances while 77 describes how 
that total abundance is distributed across the different 
modes. Indeed, we see from Eq. (|63|) that the lightest of 
the KK modes will always carry the greatest abundance. 
It then follows from its definition in Eq. ([5]) that rj indi- 
cates what fraction of the total abundance is carried by 
the excited states in the KK tower. For this reason we 
shall occasionally refer to 77 as the "tower fraction" . 

Let us first consider f2tot- At first glance, it might 
seem algebraically cumbersome to tally these individual 
mode abundances since they each have a different A- 
dependence given in Eq. (|61[) and we would need to sum 
over all of the eigenvalue solutions to the transcendental 
equation in Eq. (1551) . However, it turns out that the 
resulting spectrum of A-eigenvalues satisfies two critical 
identities [ill ]: 

J2Al = l, ]Ta^ 2 a = 1. (66) 

A A 

These identities ultimately stem from the unitary nature 
of the mapping between the KK-momentum basis and 
the mass-eigenstate basis. Thus, unless our KK tower 
experiences a staggered turn-on during the radiation- 
dominated era, summing over the individual abundances 
in Eq. (1611) is particularly simple. 

What is truly remarkable about the identities in 
Eq. (|B1))) is that they hold for all values of y. They 
are thus independent of the degree of non-diagonality 
exhibited by the KK mass matrix, and independent of 
the degree to which the corresponding KK modes are 
mixed. As a result, if we assume that f$, to, and m 
are all y- independent, it then follows that f2 t ot will be 
^-independent as well! Although dialing the value of y 
(i.e., adjusting the radius of the extra spacetime dimen- 
sion) might change the distribution of abundances across 
the infinite tower of KK states, the total abundance re- 
mains essentially fixed. 

We stress that this result holds only for those cases in 
which the KK abundances are established either instan- 
taneously, or through a staggered turn-on which takes 
place during the reheating or matter-dominated eras. By 
contrast, if these KK abundances are established through 
a staggered turn-on during the radiation-dominated era, 
the total abundance fitot today will be proportional to 
X)a ^ 1 ^ 2 A\- This quantity is not y-independent, but 
rather has the y-dependence shown in Fig. 2) As a re- 
sult, our preferred choices for parametric quantities such 
a s f<f>, to, and to would need to be altered in order to 
compensate for this effect. 

The results in Eq. (|BT|) also allow us to calculate the 
tower fraction 77 prior to the onset of significant KK 
decays. Indeed, the three different patterns for Q\ in 
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FIG. 4: The quantity A 1</2 y4 2 , plotted as a function of y. 
Note that the total dark-matter abundance Qtot prior to the 
onset of significant KK decays is proportional to this quan- 
tity if the individual KK abundances are initially established 
through a staggered turn-on during the radiation-dominated 
era. As a result, this curve also illustrates the {/-dependence 
of Qtot in this case. By contrast, for all other cases, the total 
abundance fi t ot is y-independent. 



Eq. (I61[) imply three different distributions for the abun- 
dances across the different mass eigenstates in the KK 
tower. Therefore, if we additionally assume that all of 
the states in a given tower simultaneously fall into one of 
these three cases, there will be three corresponding possi- 
ble behaviors for the "tower fraction" 77 defined in Eq. |[5} , 
viewed as a function of the non-diagonality parameter y. 
These results are shown in Fig. [5j and we see that 77 in- 
deed spans a range of 0(1) values, as desired. These re- 
sults are also highly y-dependent, illustrating that while 
adjusting y changes the total abundance only in certain 
restricted circumstances, it changes the distribution of 
these abundances quite substantially in all cases. 

It is easy to understand the overall features exhibited 
in Fig. [5J As y — > 00, we enter the four-dimensional 
limit in which virtually no abundance is carried by the 
excited KK modes. As a result, these modes become 
completely irrelevant to the dark-matter problem, and 
77 — > 0. By contrast, as y — > 0, our KK states experience 
maximal mixing, as a result of which the corresponding 
value of 77 is maximized. Since A„ ~ (n + \)/R for all 
71 in this limit, we can easily calculate these maximum 
values of 77, obtaining the result that 77 max = 1 in the case 
of instantaneous turn-on, while for the case of staggered 
turn-on during a reheating or matter-dominated era we 
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FIG. 5: The tower fraction rj after all dark-matter modes 
have "turned on" and entered the present matter-dominated 
epoch, plotted as a function of y for three different regimes of 
misalignment production: (a) instantaneous turn-on, in which 
case J1a ~ -^ 2 ^4a> 3 ) staggered turn-on during a radiation- 
dominated era, in which case Q\ ~ A 1 ' 2 ^ 2 ; and (c) staggered 
turn-on during a reheating or matter-dominated era, in which 
case fl\ ~ A 2 . In each case we see that 77 — > as y — > 00, 
while rj approaches a fixed maximum value ?7 max as y — > 0. 



have 



1-4 



1 



E 



^(» + l/2)> 
« 0.189 , 



(66) 



and for the case of staggered turn-on during a radiation- 
dominated era we have 



= 1 - 2\/2 



n=0 v 

2\/2 



1 



(2V2-l)C(3/2) 



n+ 1/2) 3 / 2 

0.408 (67) 



where £ denotes the Riemann zeta-function. All three of 
these limiting values are evident in Fig. [5j 



2. Decay widths Y\ 



M = and a non- vanishing brane mass m 7^ 0. As 
we have seen, this has proven sufficient to allow us to 
determine not only the Kaluza-Klein "spectroscopy" of 
our dark towers but also the corresponding cosmological 
mode abundances that emerge from misalignment pro- 
duction. In some sense, it is remarkable that these results 
rely on such minimal assumptions; indeed, this happy 
fact explains why our results thus far are extremely gen- 
eral and can be expected to hold for all bulk fields for 
which M = and m ^ 0. 

However, in order to discuss the decay widths of these 
KK modes into Standard-Model states, we shall require 
further information concerning the couplings between the 
five-dimensional bulk field $ and the four-dimensional 
Standard-Model states. In other words, we shall require 
further information concerning the interaction terms that 
might appear within C- m t in Eq. (|31[) and thereby be- 
come part of our four-dimensional effective Lagrangian. 
Of course, it is most dangerous for the consistency and 
phenomenological viability of our dynamical dark-matter 
scenario if these decay widths are too large. Our conser- 
vative approach to this problem will therefore be to con- 
sider the worst possible scenario and determine how large 
these decay widths might be. Since the largest decay 
widths will generally arise from the operators of lowest 
possible dimension within C- m t , the first step in our anal- 
ysis is therefore to determine what forms such operators 
might take. 

Imposing Lorentz invariance and invariance under all 
Standard-Model gauge symmetries, we find that our op- 
tions are fairly limited. The separate brane/bulk struc- 
ture of this dark-matter setup requires that our bulk 
dark-matter field $ be a singlet under all Standard- Model 
gauge symmetries. This implies that any combinations of 
brane fields to which $ couples must be gauge invariant 
by themselves. Moreover, in order to restrict ourselves 
to operators of lowest possible dimensionality, we shall 
consider operators which are at most linear in $. We 
shall also assume for this discussion that $ is real. Let- 
ting ip denote a generic Standard- Model fcrmion and F^ 
denote a generic Standard-Model field strength for any 
gauge group, we then find that the operators of lowest 
possible dimensionality come in two groups. If $ is CP- 
even, the lowest-dimension operators which may appear 
in £j nt take the form 



1 



,3/2 
./ $ 



1 



,3/2 
,/<I> 



(68) 



where /$ is the five-dimensional mass scale associated 
with $ which originally appeared in Eq. (|47j) . By con- 
trast, if $ is CP-odd, the lowest-dimension operators 
which may appear in £j nt take the form 



Next, we turn to the decay widths T\ which can be 
expected in such a scenario. 

Up to this point, we have assumed nothing more than 
that the bulk field in our setup has a vanishing bulk mass 
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f • 
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(69) 



where F^ u ~ e^ vpa F pa . These groups of operators then 
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respectively give rise to the four-dimensional couplings 



1 



(70) 



and 



X^O^WV, \<t>'F^ (71) 

where <j>', as defined in Eq. (|45|) . is the projection of $ 
onto the Standard-Model brane, and where is defined 
in Eq. (Pf. 

This list exhausts the possible dimension-five opera- 
tors. It is encouraging that we see among this list of pos- 
sible operators the standard moduli and axion couplings 
— indeed, in the CP-even case we can even regard the lin- 
ear prefactor <f)' / as the leading term of an exponential 
prefactor exp(</>'/ /^), and thereby recognize the standard 
dilaton coupling in string theory. Thus, this list of opera- 
tors includes most of our cases of phenomenological inter- 
est. At first glance, it might seem that operators of even 
lower dimension could be constructed — e.g., Qipip and 
^ipj 5 ip. However, such operators are not gauge invari- 
ant because all of the fermions ip in the Standard Model 
are chiral. Likewise, dimension-four operators such as 
|$| 2 |i/ 1 2 are also forbidden as they would violate the shift 
symmetry under which $ — > $ + constant. 

Although kinematic effects favor the di-photon de- 
cay mode (j)' — > 77, the four-dimensional couplings in 
Eqs. (1701) and (|TTj) all lead to SM decay rates of the same 
parametric order. Standard calculations then lead to an 
overall decay width 



A 



A 3 

f(j> 



(72) 



where we have substituted Eq. (14T>1) in the final step. Use 
of Eqs. (|62|) then leads to the large- A behavior T\ ~ A 3 
as well as the small- A behavior ~ A 5 . 

Before concluding our discussion of the decay widths, 
it is important to note that there will generally exist 
many competing decay modes for our KK states which 
do not exclusively involve Standard-Model particles as 
end-products. One example includes intra-ensemble de- 
cays {i.e., decays within the KK tower, from heavier KK 
states to lighter KK states); indeed, this possibility will 
be discussed in general terms in the Appendix. Also, in 
cases involving multiple fields in the bulk, it is possible 
for bulk KK states of one species to decay to bulk states 
of another species. 

While such decays can be important on a number of 
cosmological and phenomenological levels, they gener- 
ally do not significantly diminish the abundance of what 
might be termed "dark matter" or increase the corre- 
sponding abundance of what might be termed "visible 
matter" . Moreover, it is often the case that such decays 
are significantly suppressed relative to the KK decays 



that proceed directly to Standard-Model brane states. 
Such suppressions can occur for a variety of reasons, 
some of which depend on the fact that physics in the 
bulk is often governed directly by (and therefore sup- 
pressed by) the gravitational Planck scale, and some of 
which are the consequences of extra restrictive symme- 
tries which exist purely in the bulk and which therefore 
do not apply to decays of bulk fields into brane fields. 
Of course, a detailed analysis of this question requires 
specifying a particular bulk field, along with a complete 
Lagrangian for the theory including its gravitational in- 
teractions. While such an analysis is beyond the scope 
of this theoretical overview, an analysis of this sort does 
appear in Rcfs. @, 0] where it is shown that such decays 
are indeed greatly suppressed in a specific realistic model 
of dynamical dark matter. This result therefore confirms 
our general expectations in one specific example. 

We have therefore assumed in this paper that the pri- 
mary decay mode for each KK bulk mode is directly into 
a Standard-Model brane state. However, as discussed in 
the Appendix, our dynamical dark-matter scenario can 
easily be generalized to accommodate more complex de- 
cay channels if this should ultimately prove appropriate 
in a given situation. 



3. Balancing lifetimes against abundances 

Having calculated the spectrum of abundances fl\ and 
the spectrum of decay widths T\ across our KK tower, we 
can now see exactly how KK towers manage to balance 
lifetimes against abundances. Combining the results in 
Eqs. (|63|) and (|72|) . we find that for large A our KK tow- 
ers indeed always obey a balancing equation of the form 
anticipated in Eq. 



instantaneous 
staggered (RD era) 
staggered (RH/MD eras) 



.2/3 
A 

J7>r A ' " ~ constant 

A 



^aTa' " ~ constant 

.7/6 



fl\T 4 J 3 ~ constant 



(73) 



In deed, thi s asymptotic behavior holds for A 
y/l + ir 2 /y 2 . Thus, we see that our KK towers succeed in 
balancing lifetimes against abundances in a robust fash- 
ion, regardless of the particular type of "turn-on" they 
experience and regardless of the cosmological era during 
which this "turn-on" takes place. 

While Eq. (f73|) describes the crucial asymptotic behav- 
ior at the "top" of the KK tower, a similar set of relations 
describes the "bottom" of each KK tower. Combining 
the small-A behavior in Eq. (|64p with the small-A result 
Ta ~ A 5 leads to the relations 



instantaneous 
staggered (RD era) 
staggered (RH/MD eras) 



constant 



.3/10 
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i7Ar A /iu ~ constant 

.2/5 



(74) 
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for A <C y/T + ir 2 /y 2 . We stress, however, that this be- 
havior is relevant only at the bottom of a KK tower, and 
only for relatively small y. Indeed, regardless of the value 
of y, the behavior of the abundances and decay widths al- 
ways eventually shifts to satisfy the relations in Eq. (|T3[) 
as we pass to higher and higher modes in a given KK 
tower. 

Given these results for the abundances fl\ and de- 
cay widths T\, we can now calculate the general (a, (3) 
scaling coefficients that appear in Eq. (fT5|) . These re- 
sults also enable us to deduce an "effective" equation 
of state for our ensemble of decaying dark-matter KK 
components. The values of a, of course, are directly ev- 
ident from Eq. ([73)) for large A and from Eq. (|74|) for 
small A. Likewise, since the states in our KK tower are 
nearly equally spaced throughout the tower, we know 
that the density of states per unit A is essentially A- 
independent: n\ ~ A . Per unit of F, this translates 
into n r ~ nxldT/dX]- 1 ~ T^~ x ^ x for T ~ X x . We thus 
have (3 = —2/3 for large A, and (3 = —4/5 for small A. 

We therefore conclude that for large A, a general KK 
tower has the scaling coefficients 

{(—2/3,-2/3) instantaneous 
(-7/6, -2/3) staggered (RD era) 
(-4/3, -2/3) staggered (RH/MD eras) . 

(75) 

By contrast, for small A, these results are modified to 
become 

{(0,-4/5) instantaneous 
(-3/10, -4/5) staggered (RD era) 
(-2/5, -4/5) staggered (RH/MD eras) . 

(76) 

Given these (a, j3) scaling coefficients, we can also 
calculate the effective equation-of-state function w e g(t) 
which describes the collective effects of the decays of the 
individual modes along the KK tower. Indeed, as we 
have seen in Sect. II, the behavior of this function w e g (t) 
depends critically on the value of the sum x = a + ft. 
However, given the results in Eqs. (|75j) and (|76|) . we can 
easily tabulate the values of x for the different cases under 
study, obtaining the results shown in Table fl] As we see 
from Table U most of the x- values for a general KK tower 
tend to cluster near x < — 1. This is remarkable, since we 
have already shown in Sect. II that this is precisely the 
range for x which is preferred phenomenologically . We 
thus see that a KK tower indeed serves as an excellent 
realization of dynamical dark matter. 

One feature which emerges from Table Q] is that re- 
gardless of the turn-on behavior of the individual modes, 
the value of x generally decreases as we pass from the 
large-A regime to the small-A regime. This generally cor- 
responds to passing from early times (during which the 
decays of the heavier KK modes dominate the physics) to 
later times (during which only the lighter KK modes are 
still present). Indeed, this transition typically occurs for 
values of A ~ y/l + ir 2 /y 2 which decrease as a function 
of y. Thus, we do not expect to see the small-A behavior 





large A 


small A 


instantaneous 
staggered (RD era) 
staggered (RH/MD eras) 


-4/3 
-11/6 

-2 


-4/5 
-11/10 
-6/5 



TABLE I: Values of the equation-of-state parameter x = a+/3 
for different portions of a general KK tower with different 
"turn-on" phenomenologies. We observe that KK towers nat- 
urally give rise to values x ^ — 1, which is precisely the range 
favored phenomenologically. 



emerge strongly except for later times in small-y scenar- 
ios. 

Needless to say, all of the above conclusions are pred- 
icated on approximations which model the KK tower 
according to certain power-law scaling behaviors. It is 
therefore natural to wonder how robust these conclusions 
actually are when compared with the results of a com- 
plete numerical calculation which uses the exact numeri- 
cal values for the eigenvalues A across the entire KK tower 
and which avoids any approximations for the coefficients 
A\ which appear in the KK mode abundances and de- 
cay widths. However, it is straightforward to perform 
such a calculation. In Fig. [BJ we plot a rescaled version 
of the total dark-matter abundance f2tot as a function 
of time during its final decay-dominated period, assum- 
ing (as in Sect. II) that these decays occur during the 
present matter-dominated cosmological era. Each panel 
in Fig. |5] corresponds to one of the three different cases 
that describe how the individual abundances in the KK 
tower might have been established; indeed, following the 
results in Eq. (joTj) , this "rescaled" fitot is defined in each 
case as J2\ ^ A\X\ with k = 2 (first panel of Fig. [5]), 
k = 1/2 (second panel), and k — (third panel). More- 
over, in making these plots, we have assumed that each 
KK state decays instantaneously at t = t\ = T^ 1 so that 
the contributions from individual states will be readily 
discernible. This is tantamount to approximating X\ in 
Eq. (J55j) as X x (t) « QQT^ 1 -*). Finally, in order to com- 
pare curves with different values of y, an overall normal- 
ization for the time axis for each curve has been chosen 
such that the time t is expressed in units of , where 
L is the decay width of the lightest KK mass eigenstate. 
As a result the curves in Fig. [H] share a common location 
at which f2 t ot ultimately vanishes in each case, signifying 
the eventual decay of the final, lightest state in the KK 
tower. 

Although the total dark-matter abundance f^tot ulti- 
mately vanishes at log(Loi) = for all curves, we see 
that the overall time-evolution of Qtot as we approach 
this vanishing point is highly y-dependent. As y — > oo, 
we see that f2 to t remains constant until it experiences 
a single, sudden, complete decay; this of course corre- 
sponds to the traditional scenario of a single dark-matter 
particle. By contrast, for smaller values of y, we see that 
multiple modes with different decay widths carry the to- 
tal dark-matter abundance fltot', as a result, the resulting 
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FIG. 6: The (rescaled) total dark-matter abundance fitot, plotted as a function of time for each of the three cases relevant for 
a general KK tower. In each panel, the uppermost curve corresponds to y = 10 and the successively lower curves correspond 
to y = 3, y = 1, y = 0.5, and y = 0.1. In order to compare curves with different values of y, we have plotted log 10 (Toi) on the 
horizontal time axis, where To is the decay width of the lightest KK mass eigenstate associated with each curve. This ensures 
that the curves share a common horizontal location at which Q to t vanishes in each case, signifying the decay of the final, lightest 
state in the KK tower. We see from these results that overall shape of the time-dependence of 57 to t is highly {/-dependent: the 
y — > oo limit corresponds to the usual scenario of a single dark-matter particle (with Sl to t remaining essentially constant until 
this particle decays), while smaller values of y correspond to situations in which Q to t is distributed across multiple KK modes 
with different decay widths. It is this property which leads to a time-dependent S7tot and thus a non-trivial "effective" equation 
of state for the dark KK tower. 




staggered (RD era) 




staggered (RH/MD eras) 




FIG. 7: The tower fraction r), plotted as a function of time for each of the three cases relevant for a general KK tower. In each 
panel, the lowest curve corresponds to y = 10 (not visible in the center and right panels) and the successively higher curves 
correspond to y = 3, y = 1, y = 0.5, and y — 0.1. Thus, these values for n directly correspond to the values of fitot plotted in 
Fig. [6] In general, we observe that r\ increases with decreasing y, and that the maximum value of r\ shown for each curve is 
consistent with the results of Fig. [5] 



transition of fitot from its maximum value to zero is more 
gentle. In all cases the quantity 1 — rj indicates the rela- 
tive size of this final "last step" down to fitot = 0; note 
that the results for rj implicitly shown here in terms of 
the relative final step size are consistent with those shown 
in Fig. [5] Likewise, the initial values of fitot also confirm 
our expectations discussed earlier: the instantaneous and 
staggered (RH/MD) cases have initial values at f2 to t = 1, 
in accordance with Eq. (|65|) for all y, while the initial 
values shown on the second panel are y-dependent and 
correspond to the values shown in Fig. |4j 

Using the results for fl t ot(t) shown in Fig. [BJ we can 
now proceed to calculate the corresponding tower frac- 



tions rj{t). The results are shown in Fig. [7J As expected, 
we see that rj increases with decreasing y. Moreover, we 
can now see directly that the maximum value of rj shown 
for each curve is consistent with the results of Fig. [S] 

Using the results for Q to t(t) shown in Fig. [5J we can 
also proceed to calculate the corresponding equation-of- 
state function w e ff(t) which follows from the definition 
in Eq. ([7]). Equivalently, we can use the results for r](t) 
shown in Fig. [JJ along with the definition in Eq. ©. In 
either case, the results are shown in Fig. [HI In passing 
from Figs. [S] and [JJ to Fig. [SJ we have calculated loga- 
rithmic slopes numerically for each successive KK decay 
event and then plotted a continuous function. It is clear 
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FIG. 8: The effective dark-matter equation-of-state parameter w c s, piotted as a function of time for each of the three cases 
relevant for a general KK tower. In each panel, the lowest curve corresponds to y = 10 and the successively higher curves 
correspond to y = 3, y = 1, y = 0.5, and y — 0.1. Thus, these values for w e B directly correspond to the values of f2 to t plotted 
in Fig. [6l or equivalently the values of r\ plotted in Fig. [7] In all cases, we see that w e s — > as t — > 0. Note that although the 
values of f2 to t were plotted in Fig. [6] only up to an overall rescaling factor, the values of w c k plotted here are insensitive to this 
rescaling and are thus meaningful on an absolute scale. We thus see that w c g never exceeds 0.1 for a general KK tower, and is 
generally much smaller. 



that the results in Fig. [5] are in complete agreement with 
our general expectations for w e g from Sect. II: in each 
case we observe the general tendency that w c g — > as 
t — > 0, and likewise in all but one case w c g approaches 
a pole at t = Tq 1 (corresponding to the fact that fitot 
hits zero upon the decay of the final, lightest dark-matter 
mode in the KK tower). 

These results are also in agreement with our expecta- 
tions based on the KK scaling coefficients in Table HI As y 
decreases, we see the emergence of a definite shift in the 
behavior of w e {f(t) as we transition from early times to 
later times. This corresponds to the shift from large-A be- 
havior to small- A behavior in Table fl] Indeed, in the case 
of an "instantaneous" turn-on, we even observe that our 
function w c g{t) develops a slight maximum for smaller 
values of y, shifting from increasing behavior to slightly 
decreasing behavior. This change in the slope of w c s(t) 
for this particular situation is directly correlated with the 
fact that the value of x in Table Q] shifts from x < — 1 to 
x > — 1. We see from both Table U and Fig. [5] that this is 
the only case in which such interesting behavior occurs. 

At first glance, it might seem surprising that we are 
able to obtain effective equation-of-state functions w c f{(t) 
which depend on x, but which are otherwise universal 
when plotted versus Tot. Indeed, w Q g(t) depends on a 
number of parameters: not just the dimensionless expo- 
nents a and /3 in Eq. (|15p . but also dimensionful quanti- 
ties such as the leading coefficients A and B in Eq. (fT5|) 
as well as physical parameters such as JIcdm which are 
involved in setting a boundary condition for tttot- Indeed, 
all of these parameters appear in our approximate results 
for fitot(i) in Eqs. (fT8|) and (fT9|) . However, the important 
point is that while f2tot depends on all of these dimension- 
ful quantities somewhat independently, w e g(t) depends 
on them in only one particular combination. This was 



already apparent in Eqs. (|2U)) and (|2U)) . where the com- 
bination in question was nothing but w* = w c s(t now ). 

Of course, the results for w e s(t) in Eqs. (|2"D|) and (j2"5)l 
were respectively derived from the results for O to t (t) in 
Eqs. ([IH]) and ([TO)) , and these in turn were realized by 
taking our boundary condition to be fltot(t now ) = f^cDM- 
However, we can equivalently write our boundary condi- 
tion in the form f2 t ot(l/ro) = 0, where we are assum- 
ing that each KK mode with mass A decays promptly at 
t = t\ = r^ 1 , with Tq denoting the width of the lightest 
mass eigenmode. Following the same algebraic manip- 
ulations as in Sect. II then leads to equations of state 
which are written in terms of Tq rather than : 



[ i(x + l)[l-(r ir +1 ] 1 forx^-1 
w e s(t) = < 

[ i(roi)- 1 for x = -1 

(77) 



where x = a + f3. Thus, when expressed in terms of the 
dimensionless time variable Tot as in Fig. |8l our w c s(t)- 
functions are indeed universal, depending only on x. 

It is important to bear in mind that in this section we 
have made only minimal assumptions concerning the pre- 
cise nature of this KK tower or the identity of the fields 
it represents. We therefore expect that all of the features 
we have discussed in this section will hold quite gener- 
ally, regardless of the identity of the particular field(s) 
which happen to populate the bulk and constitute the 
KK tower. 
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IV. DYNAMICAL DARK MATTER: NOVEL 
SIGNATURES AND PHENOMENOLOGICAL 
CONSTRAINTS 

Having described the general theoretical structure of 
our dynamical dark-matter framework, we now present 
several additional features of this framework which are 
likely to be of importance in enabling this framework to 
satisfy phenomenological constraints. As in other sec- 
tions, our discussion here will be restricted to broad, 
model-independent themes, and we shall present a de- 
tailed phenomenological analysis of one specific dynami- 
cal dark- matter scenario in Refs. 0,0. 

It turns out that there are three phenomenological fea- 
tures which are unique to dynamical dark matter, and 
which under certain circumstances might be taken as sig- 
natures (or even "smoking guns") of the overall frame- 
work. 

• No well-defined dark-matter mass or cross-section: 
First, since the dark-matter "candidate" within the 
dynamical dark-matter framework is not a single 
particle, but rather an ensemble of particles, it 
does not have a specific mass or cross-section asso- 
ciated with it. This represents a marked difference 
relative to most other dark-matter proposals, and 
implies that the kinematics associated with the 
production and decay of dynamical dark matter 
is likely to be quite different from that associated 
with more traditional single-component dark 
matter. This could have dramatic consequences 
for collider phenomenology (and potentially for 
direct detection), and may also have a number of 
cosmological implications. 

• Coupling suppression for light modes: Second, it 
is almost inevitable that the eventual phenomeno- 
logical success or failure of specific dynamical dark- 
matter scenarios involving large extra dimensions 
will ultimately rest in part on the couplings be- 
tween the dynamical dark matter in the bulk and 
the Standard-Model states on the brane. Such cou- 
plings are of critical importance because they gov- 
ern the degree to which this dark matter might 
be "visible" to the Standard Model. From intu- 
itions based on studies of KK graviton dynamics, 
one might suspect that all mass eigenstates in the 
bulk would couple to the brane with equal strength. 
However, it turns out that the opposite is true for 
theories with a brane mass: while the most of the 
couplings between the fields on the brane and the 
mass eigenstates in the bulk are indeed uniform, 
the couplings between the brane and the lightest 
mass eigenstates in the bulk are significantly sup- 
pressed. We have already seen this behavior in 
Fig. [3l where we plot the coupling matrix element 
{4>\W) = ^ 2 ^I as a function of A: although this 
coupling always reaches an asymptotic value for 
sufficiently large A, this coupling is significantly 



suppressed for small A. Indeed, the magnitude of 
this suppression can be controlled by varying the 
non-diagonality parameter y. 

This suppression feature is of immense phenomeno- 
logical importance, since the couplings of the light- 
est dark-matter eigenstates to the brane are pre- 
cisely those which are the most dangerous for 
the viability of our dynamical dark-matter frame- 
work. Thus, this suppression feature can be very 
important in relaxing numerous phenomenological 
bounds on dynamical dark matter, and thereby 
constitutes an unexpected additional effect which 
can help dynamical dark matter stay dark despite 
its multitude of states. 

• Decoherence: Finally, there is an additional fea- 
ture associated with dynamical dark matter in large 
extra dimensions which can play an important role 
in its phenomenology: this is the phenomenon of 
"decoherence" As we have seen in Sect. Ill, 

only one particular linear combination of bulk dark- 
matter modes 4> n can couple to the brane: this 
is the linear combination cf>'. However, once 0' is 
created through an interaction with the brane, it 
rapidly decoheres as it propagates because it is not 
a mass eigenstate. 

One way to understand this decoherence involves 
simple quantum mechanics: because <p' consists 
of a huge number of different mass eigenstates, 
and because the masses of these eigenstates are 
generally not related to each other through ratio- 
nal multiplicative factors, the different mass eigen- 
states fall out of phase with each other under time- 
evolution and will not reconstitute <f>' within fi- 
nite time. Thus, they cannot couple to the brane 
at later times, and essentially become "invisible" 
as far as physics on the brane is concerned. An- 
other (quantum field-theoretic) way to describe the 
same phenomenon is simply that the amplitude as- 
sociated with any process that involves the pro- 
duction and subsequent detection of dark matter 
on the brane will have multiple individual contri- 
butions, each associated with the propagation of 
an intermediate state consisting of an individual 
dark-matter component. However, because these 
individual components have different masses, their 
corresponding amplitudes accrue different phases. 
These amplitudes therefore destructively interfere 
within the calculation of any cross-section sum. 

This decoherence phenomenon can have impor- 
tant phenomenological consequences. Indeed, de- 
coherence generically induces a suppression of the 
cross-section for any process involving virtual dark- 
matter particles by a factor of N, where N is the 
number of such particles being exchanged. This, 
then, is yet another mechanism which helps dy- 
namical dark matter stay dark. We emphasize that 
this feature is not specifically extra-dimensional; 
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it applies to any dark-matter framework in which 
the dark matter has many components of different 
masses, and in which only a specific linear combina- 
tion of those components can couple to Standard- 
Model states. 

Needless to say, dynamical dark matter must ulti- 
mately be subjected to all of the phenomenological 
bounds and constraints that apply to more traditional 
dark-matter candidates. However, because dynamical 
dark matter consists of a vast ensemble of individual 
states which are not necessarily stable on cosmological 
time scales, many of these constraints take unusual forms 
in this context. We shall therefore now provide a quick 
overview of the different classes of laboratory, astrophys- 
ical, and cosmological constraints which apply to dark 
matter in general, and then indicate the forms they can 
be expected to take within the context of dynamical 
dark matter. Once again, we emphasize that our goal 
here is merely to provide a model-independent theoret- 
ical overview in which we restrict ourselves to address- 
ing a single question: for each class of constraints that 
apply to theories of dark matter, what combinations of 
parameters are bounded in the traditional framework and 
how do these combinations translate into our dynami- 
cal dark-matter framework? Explicit details concerning 
a particular dynamical dark-matter scenario can be found 
in Refs. 

Broadly speaking, there are four classes of constraints 
which apply to any candidate theory of dark matter. 

• First, there are general constraints on the relic 
dark-matter abundance and on the dark-matter 
equation of state. The constraints on the dark- 
matter abundance f2tot(i) are similar to those 
which apply in traditional dark-matter scenarios: 
^tot(^now) must match the observed dark-matter 
relic abundance Hcdm; the dark-matter ensemble 
must not cause the universe to become matter- 
dominated too early; etc. However, our sce- 
nario differs from traditional dark-matter scenar- 
ios in that it gives rise to an equation-of state- 
parameter w c g which can be different from zero and 
which generally exhibits a non-trivial time depen- 
dence. Astrophysical and cosmological considera- 
tions therefore imply additional constraints on our 
equation-of-state function w c s(t), or equivalently 
on the scaling coefficients (a, (3). 

• Second, there are constraints on dark matter which 
derive from physical processes in which dark matter 
is produced through its interactions with Standard- 
Model particles but not subsequently detected. For 
example, there are collider constraints on processes 
in which dark-matter particles manifest themselves 
as missing energy. Furthermore, if the dark-matter 
candidates in question are sufficiently light, addi- 
tional constraints can be derived from limits on 
dark-matter production by astrophysical sources. 



For example, dark-matter particles produced in 
stars and supernovae can carry away energy from 
these sources very efficiently. This can lead to ob- 
servable effects on stellar lifetimes, energy-loss rates 
from supernovae, etc. Indeed, observational lim- 
its on the magnitudes of these effects imply strin- 
gent bounds on any light particle whose interac- 
tions with the Standard-Model fields are highly 
suppressed. 

To see how such considerations constrain the pa- 
rameters of a generic dark-matter model, let us 
consider a traditional single-component scenario in 
which the dark-matter candidate resides in a hid- 
den sector. The dominant interaction between the 
dark sector and the Standard Model in such sce- 
narios occurs through non-renormalizable opera- 
tors O n of mass dimension n, suppressed by inverse 
powers of some large mass scale A associated with 
the cutoff of the theory. For example, A might be 
an effective Planck scale Mp in the case of dark- 
matter candidates associated with gravity, or a par- 
ticle mass Mr in the case of candidates such as a 
right-handed neutrino, or a dynamical scale such 
as the Peccei-Quinn scale /pq in the case of ax- 
ions. The cross-section for dark-matter production 
will therefore be suppressed by a factor of A 2 ( 4 ~ n \ 
and constraints in this class thus ultimately become 
constraints on A. 

For example, in the specific dynamical dark-matter 
scenario presented in Sect. Ill, the leading oper- 
ators have mass-dimension five, and A is equated 
with the suppression scale f^. Thus, constraints in 
this class ultimately become bounds on Or, 
phrased directly in terms of the decay widths and 
abundances which are the bedrock of our scenario, 
these constraints yield bounds on A~ 3 Fa, where 
the sum over mass eigenstates includes only those 
states which are kinematically relevant for the pro- 
cess in question. 

We conclude, then, that constraints of this type 
tend to place bounds on the particular combina- 
tion ~ 3 Ta- However, this quantity is signif- 
icantly affected by the coupling-suppression effect 
discussed above. As a result, such bounds can often 
turn out to be considerably weaker that one might 
imagine at first glance. 

• Third, there are constraints that arise from situa- 
tions in which dark matter is produced through its 
interactions with Standard-Model particles and is 
then subsequently detected (either directly or indi- 
rectly) via those same interactions. Here we have in 
mind not only astrophysical production with sub- 
sequent detection on earth, but also any process in- 
volving virtual dark-matter particles. Which phys- 
ical processes of this sort serve to constrain a par- 
ticular dark-matter particle are extremely model- 
specific. Axions and other similar particles, for 
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example, are constrained primarily by helioscope 
searches, microwave-cavity experiments, etc.; other 
particles are more stringently constrained by col- 
lider limits, and so forth. Nevertheless, a few 
generic observations can be made. 

If we assume, as above, that the dark matter re- 
sides in a hidden sector, it then follows that the 
cross-sections for processes of this sort are propor- 
tional to A 4 ( 4 -™). Thus, once again, limits on such 
cross-sections ultimately become bounds on A. For 
example, in the specific dynamical dark-matter sce- 
nario presented in Sect. Ill, they become bounds on 
or equivalently on the quantity (J2\ ^~ 3 F\) 2 . 
In terms of overall mass scales, we might approx- 
imate this quantity as X)a^ _6 ^A' but we must 
also bear in mind that the cross-terms within such 
products can be significant. Indeed, these are pre- 
cisely the situations in which the decoherence phe- 
nomenon discussed above can play a role. Thus, 
these constraints might also turn out to be signifi- 
cantly weaker than they might at first sight appear. 

• Finally, there are constraints on dark-matter decays 
and annihilations. As far as decays are concerned, 
we have in mind constraints such as those from big- 
bang nucleosynthesis, measurements of the cosmic 
microwave background, observations of the diffuse 
X-ray and gamma-ray backgrounds, etc. The ba- 
sic idea behind all of these constraints is that the 
decays of a cosmological population of dark-matter 
particles can result in measurable deviations from 
the standard cosmology at times t > 1 s, or leave 
(unobserved) imprints on these backgrounds. In 
situations in which dark-matter annihilation cross- 
sections are sufficiently large, various additional 
limits (such as those from typical indirect-detection 
methods) would also apply. 

Let us focus on those constraints related to dark- 
matter decays, as these are generic to dynamical 
dark-matter scenarios. (By contrast, constraints 
related to dark-matter annihilation tend to be 
somewhat model-dependent, and indeed often do 
not apply.) Roughly speaking, in a traditional 
single-component dark-matter scenario, such dark- 
matter decay constraints tend to place bounds on 
the product fl x T x , where Q x and T x are the abun- 
dance and decay width of our dark-matter field Xi 
suitably evaluated during the appropriate cosmo- 
logical period. In a dynamical dark-matter frame- 
work, by contrast, this now becomes a constraint 
on 

(r) ee £n A r A (78) 

A 

where again our abundances and widths are to 
be evaluated during the cosmological epoch dur- 
ing which decays can contribute to the effect in 
question (the disruption of BBN, distortions of the 



CMB, etc.). It is important to note that this de- 
pendence on time effectively truncates the range of 
the above sum to those states whose lifetimes fall 
roughly within the time scale associated with that 
epoch. Put another way, only those states whose 
masses lie within certain characteristic ranges con- 
tribute to the sum. Of course, at a mathematical 
level, the behavior of this sum ultimately depends 
critically on the balancing relations that happen to 
hold across the entire dark-matter ensemble. Note 
that these arguments will be addressed more rigor- 
ously in Refs. 

Finally, for completeness, we also mention two further 
classes of constraints which must also be borne in mind. 
Unlike the previous constraints, these are substantially 
more model-dependent. 

• First, there can be phenomenological bounds that 
accompany (and are therefore specific to) partic- 
ular realizations of dynamical dark matter. For 
example, we have seen that an infinite tower of 
Kaluza-Klein states furnishes an excellent realiza- 
tion of dynamical dark matter. However, the 
extra-dimension brane/bulk framework brings with 
it a whole host of additional bounds and con- 
straints, some of which come from the fact that 
we are now attempting to do standard physics 
within such a context (e.g., the need for a late- 
time- reheating (LTR) cosmology [l0|), and oth- 
ers of which place bounds on the context itself 
(e.g., Eotvos-type or Cavendish-type "fifth-force" 
experiments which restrict the allowed sizes of the 
extra dimensions). Such constraints are clearly 
highly model-dependent, and frequently they are 
also wholly independent of the general dynamical 
dark-matter framework. 

• Finally, there can also be constraints that arise sim- 
ply for reasons of theoretical self-consistency. For 
example, if we assume (as we have done here) that 
our initial dark-matter abundances are determined 
through misalignment production, then we must 
insist that misalignment production indeed dom- 
inates over other production mechanisms such as 
thermal production. This, of course, yields a non- 
trivial constraint on the parameters of the model. 
Likewise, the assumption that dynamical dark mat- 
ter in the bulk decays preferentially to Standard- 
Model states on the brane, rather than to other dy- 
namical dark-matter states in the bulk, implies yet 
another self-consistency requirement. Once again, 
however, constraints in this class tend to be highly 
model-dependent and therefore do not represent 
generic constraints on the dynamical dark-matter 
framework. 

This is clearly a fairly long list of constraints, and one 
must not minimize the impact that they can have in rul- 
ing out specific dark- matter proposals. In Refs. [f| @], 
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however, we shall study one particular realization of dy- 
namical dark matter, and we shall exhaustively work 
through all of the constraints relevant for this partic- 
ular realization. We shall find that for this particular 
scenario, our dynamical dark-matter framework indeed 
survives all known laboratory, astrophysical, and cosmo- 
logical constraints. This will thereby furnish us with an 
"existence proof" that dynamical dark-matter can indeed 
be a viable, alternative framework for addressing the cen- 
tral questions in dark-matter physics. 

Finally, we remark that for some purposes it may 
also be interesting to consider the phenomenology as- 
sociated with our overall dark-matter framework when 
the individual component decay widths Tj are extremely 
small. Of course, in the actual limit I\ — > we know 
that f2tot approaches a constant in the final, matter- 
dominated era; likewise, w c ff — > 0. In this respect, this 
limit of our framework begins to resemble a traditional 
non-dynamical dark-matter framework, thereby allowing 
us to evade many of the most stringent phenomenological 
constraints coming from dark-matter decays in the early 
universe. However, even in this limit, our framework 
nevertheless continues to retain those distinctive features 
which stem from its underlying multi-component nature. 
For example, we can still have rj ^ 0. We also still have 
the possibility of staggered turn-ons as well as the possi- 
bility of coupling suppression for light modes. We even 
continue to have decoherence, even though many of the 
processes for which decoherence is most phenomenolog- 
ically relevant will already tend to be suppressed in the 
r,; — > limit. And perhaps most importantly, the dark 
matter in our framework will continue to evade simple 
characterization in terms of a single well-defined mass or 
cross-section. 

Needless to say, we are not particularly interested in 
the limit — > 0. Indeed, we regard the dynamical as- 
pects of our dark-matter framework to be among its most 
intriguing features and key signatures. However, the free- 
dom to tune the values of Ti relative to the other dimen- 
sional parameters in our framework is important from a 
theoretical standpoint because it illustrates that our over- 
all dark-matter framework possesses a means of "dialing" 
the scale associated with its dynamical aspects relative to 
those associated with its multi-component aspects. This 
is particularly relevant because the dynamical aspects of 
our framework are often ultimately subject to an entirely 
different set of phenomenological bounds and constraints 
than those governing its multi-component aspects. Thus, 
the freedom to independently adjust the scales associated 
with these different aspects of our dark-matter framework 
gives this framework an added flexibility when it comes 
to satisfying many of the phenomenological bounds dis- 
cussed above. 

Of course, within the particular higher-dimensional 
brane/bulk context discussed in Sect. Ill, this free- 
dom may initially appear to be lacking: a single five- 
dimensional mass scale /$ governs not only the magni- 
tudes of the abundances of individual dark-matter com- 



ponents but also the magnitudes of their corresponding 
decay widths. Indeed, in particular realizations of this 
framework, even the brane mass m can be tied to /$, 
and we shall see an explicit example of this in Ref. [6|. 
However, there is in principle no reason why the mass 
scale /$ which appears in Eq. (|47p and which ultimately 
sets the scale for dark-matter abundances needs to be the 
same as the mass scale /$ which appears in Eqs. (|68l) and 
(|69|) and which ultimately sets the scale for decay widths. 
Indeed, identifying these two quantities is merely a min- 
imal assumption about the energy scales in our higher- 
dimensional theory, and we can easily envision more com- 
plex scenarios in which these two mass scales are signifi- 
cantly different. 



V. CONCLUSIONS AND DISCUSSION 

In this paper, we have introduced a new framework for 
dark-matter physics which we call "dynamical dark mat- 
ter" . Unlike most approaches to the dark-matter problem 
which hypothesize the existence of a single, stable, dark- 
matter particle, our dynamical dark-matter framework 
may be characterized by the fact that the requirement of 
stability is replaced by a delicate balancing between cos- 
mological abundances and lifetimes across a vast ensem- 
ble of individual dark-matter components. This setup 
therefore collectively produces a time-varying cosmolog- 
ical dark-matter abundance, and decays of the different 
dark-matter components can occur continually through- 
out the evolution of the universe. 

Although our framework is quite general and need not 
be tied to a specific set of particles or theoretical mod- 
els, we have shown that one natural realization of this 
scenario consists of a tower of KK states correspond- 
ing to a single higher-dimensional field propagating in 
the bulk of large extra spacetime dimensions. Indeed, 
as we have shown, the states in such a "dark tower" 
naturally obey inverse "balancing" equations that relate 
their abundances and decay widths in just the right man- 
ner. Remarkably, this remains true even if the stability 
of the KK tower itself is entirely unprotected. Our dy- 
namical dark-matter scenario is therefore well-motivated 
both in field theory and string theory, and can even be 
used to constrain the cosmological viability of certain 
limits of string theory. We have also seen that within 
this context, dynamical dark matter generically gives 
rise to certain phenomena such as coupling suppression 
and decoherence which may help to explain why dark 
matter is dark and thus far unobserved. Such phenom- 
ena transcend those usually associated with traditional 
single-component dark matter, and may in some sense 
be viewed as unique signatures for a dark-matter fram- 
work such as ours which rests on the existence of a large 
multitude of individual dark-matter components. 

Needless to say, there are many possible generalizations 
of our basic dynamical dark-matter framework. Some of 
these apply to dynamical dark matter in general, while 
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others are more specific to realizations involving extra 
dimensions. For example, insofar as our general dark- 
matter ensemble is concerned, there are several natural 
extensions which can be contemplated. 

• Not all components within the ensemble need be 
scalars. Higher-spin fields may also be considered. 
We may even demand that our ensemble be su- 
pcrsymmctric, although there would be no obvious 
need for i?-parity within such supersymmetric ex- 
tensions as far as dark-matter considerations are 
concerned. 

• In this paper, we have examined the case of rela- 
tively simple dimension-five couplings between the 
components in the ensemble and the fields of the 
Standard Model. However, different scenarios may 
involve different coupling structures, and thus dif- 
ferent models will lead to their own distinctive phe- 
nomenologies. 

• Continuing along these lines, we have assumed in 
this paper that all of the components of our dy- 
namical dark-matter ensemble are neutral under 
the Standard-Model gauge symmetries. While this 
choice is particularly convenient, allowing the pos- 
sibility of specific realizations of our scenario in 
higher-dimensional brane/bulk Kaluza-Klein the- 
ories, there is nothing intrinsic to the dynamical 
dark-matter framework that requires this to be the 
case. In particular, some or all of the compo- 
nents of our dark-matter ensemble could have 0(1) 
SU{2) weak interactions with the Standard Model. 
This would, of course, undoubtedly tighten many 
of the phenomenological constraints on such sce- 
narios; likewise, scattering processes involving such 
dark-matter components are also generally likely 
to play an important role and would need to be in- 
cluded in the analysis along the lines discussed in 
the Appendix. However, as long as the lifetimes of 
the ensemble components are sufficiently balanced 
against their abundances, the basic features of our 
framework will remain intact. 

• In this paper, we have considered all decays of 
our ensemble components to be essentially instan- 
taneous. However, such decays really have an ex- 
ponential time-dependence. The fact that these de- 
cays have different widths can thus lead to further 
non-trivial effects on the time-dependence of the 
total dark-matter abundance associated with the 
ensemble. 

• The primary decay mode for a given dark-matter 
component within our ensemble need not always 
be directly into Standard-Model states. In par- 
ticular, it is also possible to consider decays from 
heavier ensemble components into lighter ensemble 
components. Note that in this sense, we are view- 
ing the ensemble as consisting of all states which 



are neutral under Standard-Model symmetries, in- 
cluding fields which reside in what might in more 
traditional contexts be considered a hidden sector. 
Such intra-ensemble decays could significantly alter 
the sorts of balancing equations which might arise 
across our dynamical ensemble, and as we shall dis- 
cuss in the Appendix, they can thereby modify the 
time-dependence associated with fitot, f]-, and w e s . 

• Misalignment production need not be the only 
mechanism through which the abundances of our 
individual components are initially established. 
Many other mechanisms (e.g., thermal production, 
decays arising from topological defects, etc.) also 
provide ways of populating the different compo- 
nents, and can likewise lead to different resulting 
phcnomcnologies for dynamical dark matter. In 
particular, it would be very interesting (and rele- 
vant for our overall framework) to see whether the 
correct sorts of balancing relations might arise for 
situations in which our different ensemble compo- 
nents are populated in the manner of a standard 
WIMP — i.e., by thermal freeze-out. 

• At many points in this paper, we have made as- 
sumptions that simplify our analysis. For example, 
in Sect. II we have assumed that t\, the time by 
which our staggered turn-on has ended, is less than 
t<i, the time at which significant dark- matter de- 
cays commence. Likewise, we have assumed for 
much of our analysis of abundances in Sect. Ill 
that a staggered turn-on, if it occurs, takes place 
entirely within a single epoch (either RH, RD, or 
MD). While such assumptions prove useful for ana- 
lyzing the effects of different features of our frame- 
work individually, there is nothing intrinsic to the 
dynamical dark-matter framework which requires 
that these features be separated in this way, and 
numerous extensions and combinations of these fea- 
tures are possible. 

• Although we have discussed several different signa- 
tures which arc unique to dynamical dark matter, 
it is likely that our discussion has only begun to 
scratch the surface. It would be interesting to in- 
vestigate what other kinds of signatures are also 
possible within this framework. 

• Finally, our discussion in this paper has assumed 
a standard FRW cosmology. However, it would be 
interesting to repeat this analysis for a ACDM cos- 
mology (and also for versions thereof with low re- 
heat temperatures, as appropriate for theories with 
large extra dimensions). Indeed, within such cos- 
mologies, quantities such as f2 to t will experience 
additional types of time-dependence beyond those 
considered here. 

Likewise, within the specific framework of large extra 
dimensions in which our dynamical ensemble of dark- 
matter components is represented by an infinite tower of 
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Kaluza-Klein states, there are also numerous generaliza- 
tions and extensions which may be contemplated. 

• We may consider situations involving multiple 
species of bulk fields. For example, the bulk can 
be a crowded place, consisting of a whole plethora 
of particles which are neutral under all Standard- 
Model gauge symmetries: these include gravitons 
and gravitini, axions and other axion-like particles, 
string-theory moduli, right-handed neutrinos, and 
so forth. From the point of view of physics on the 
brane, all of these states can be considered "dark 
matter" , and their contributions to quantities such 
as fi to t must all be considered within the overall 
dynamical dark-matter framework. 

• We also need not restrict ourselves to a single extra 
spacetime dimension. Multiple extra dimensions 
are also possible. 

• Likewise, our extra dimensions need not necessar- 
ily be flat. Warped extra dimensions will give rise 
to an entirely different KK spectroscopy, and as a 
result the phenomenology of dynamical dark mat- 
ter within such contexts is likely to be significantly 
different from what has been presented here. 

The above represent ideas for generalizing our overall 
dynamical dark-matter scenario. However, it may also be 
possible to use this kind of dynamical framework in order 
to address questions that go beyond dark matter per se. 
While some of these are relatively straightforward, others 
are indeed quite speculative. 

• One of the key features of dynamical dark mat- 
ter is that quantities such as f2 to t ar< 3 dynamical 
(i.e., time-dependent) in this framework — even 
during the current matter-dominated epoch dur- 
ing which the dark-matter abundance is normally 
thought to be roughly constant. It is therefore pos- 
sible that such a dynamic approach could serve as 
a starting point towards addressing the cosmic co- 
incidence problem. 

• Further along these lines, it might also be possi- 
ble to address the cosmological constant within a 
similar framework. For example, the energy den- 
sity associated with each scalar field 4>i prior to 
its "turn-on" behaves as dark energy rather than 
dark matter. Thus in this respect the cosmolog- 
ical constant in our framework is time-dependent 
as well, and this sort of dynamic cosmological con- 
stant might even persist into the current epoch if 
there continue to exist light scalar modes which 
have not yet turned on. In this case, the process of 
dark-matter decays would necessarily overlap with 
the process of a staggered turn-on, and one might 
be able to develop a consistent theory in which a 
vast ensemble of states gives rise to both dynamical 
dark matter and dynamical dark energy. 



• The framework of dynamical dark matter might 
also provide a new means of placing phenomenolog- 
ical bounds on string theory (and in particular, on 
candidate string models). After all, string models 
are typically rife with "bulk" fields — even if their 
extra dimensions are compactified at or near the 
traditional Planck scale. Some of these fields (such 
as moduli) are model-dependent: they depend 
on the particular kind of compactification geom- 
etry employed in the construction of the candidate 
string model and their masses depend on the par- 
ticular stabilization mechanisms, if any, which have 
been employed. By contrast, some of these fields 
are generically model-independent: these include 
all fields associated with the (super)gravity multi- 
plet, such as the graviton, dilaton, other higher- 
form fields, and their possible superpartners. In- 
deed, if the Standard Model is restricted to a stack 
of D-branes within a given string construction, the 
corresponding "bulk fields" include all string states 
which do not couple to those branes. 

While these fields are typically required for the self- 
consistency of the string, those that are massless 
require stabilization. Indeed, this is nothing but 
the standard moduli problem. However, depend- 
ing on the specific cosmological properties of these 
fields, it is also true that their abundances must 
necessarily be considered as contributing either to 
the total dark energy or the total dark matter of 
the universe — even after they are stabilized. An 
analysis of their cosmological effects is then likely to 
run along the lines we have presented here, and the 
cosmological viability of the underlying candidate 
string model thus necessarily becomes an issue to 
be studied within a dynamical dark-matter (or dy- 
namical dark-energy) framework. Indeed, in such 
cases our dark-matter ensemble could potentially 
include not only string KK modes, but also (a sub- 
set of) string oscillator modes and string winding 
modes. 

• One central feature of our dynamical dark-matter 
scenario is the phenomenon in which an ensemble of 
decaying "stuff" with one equation of state collec- 
tively simulates "stuff" with a different equation of 
state. For example, in the specific dynamical dark- 
matter scenario presented here, an ensemble of de- 
caying dark- matter states (i.e., each with w = 0) 
collectively simulates an effective equation of state 
with w c ff > 0. This notion of using a vast ensem- 
ble of states with one equation of state to simulate 
another is, we believe, worthy of exploration in its 
own right, independently of the specific uses for 
dark-matter physics that we have presented here. 

• Further along the above lines, it is natural to ask 
whether we could construct an ensemble of indi- 
vidual components of "stuff" with negative w. The 
decays of these components within the ensemble 
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would therefore act to increase the effective value 
of w, and perhaps even simulate w = 0. In other 
words, it is possible that dark matter might not 
even need to be comprised of matter! In some 
sense, this is the converse of the scenario we have 
presented here, in which individual matter com- 
ponents collectively produce a value for w e g which, 
though not too different from zero, is still non-zero. 
Indeed, both scenarios may represent equally legit- 
imate approaches to the dark-matter problem. 

• Pursuing these lines still further, one can even spec- 
ulate as to whether an ensemble of decaying dark- 
energy components (each with w = — 1) could sim- 
ulate dark matter (with w = 0). Indeed, such 
individual dark-energy components need be noth- 
ing more complicated than a set of scalar fields 
(j>i which decay (potentially into a hidden sector) 
prior to turning on. Such an approach might then 
"unify" dark energy and dark matter as simulta- 
neously stemming from a primordial ensemble of 
scalar fields. 

• Along entirely different lines, there is another phe- 
nomenon inherent in our dynamical dark-matter 
framework which is potentially interesting in its 
own right: this is the phenomenon (discussed in 
Sect. Ill) in which a KK tower appears to have 
periodic modings for its heavier modes, but anti- 
periodic modings for its lighter modes. As we have 
seen in Sect. Ill, this result emerges rather generi- 
cally, requiring only a bulk field that somehow ac- 
crues a non-zero brane mass. This phenomenon is 
extremely interesting, because one normally asso- 
ciates the modings of a given field with its bound- 
ary conditions around non-contractible loops in a 
topologically non-trivial space, or equivalently with 
the magnitudes of the fluxes which might penetrate 
those loops. This phenomenon therefore seems to 
suggest a mechanism by which such modings or 
fluxes might become effectively energy-dependent. 

• Finally, the general phenomenon of decoherence is 
interesting in its own right. This might be ex- 
tremely relevant for axion invisibility (see, e.g., the 
discussion in Ref. 

EH). 

and is also likely to be of 
more general applicability. Indeed, this might pro- 
vide an interesting approach to the moduli problem 
in string theory, and explain why moduli such as 
the dilaton are not observed. 

We see, then, that our dynamical dark-matter frame- 
work appears to be pregnant with numerous possibili- 
ties for extension and generalization. However, even as a 
framework for dark-matter physics, we caution that our 
presentation here has been limited to only the broadest 
model-independent theoretical aspects and features. In 
particular, it still remains to choose a specific realization 
of this scenario in terms of a particular species of bulk 
field, and examine the phenomenological consequences of 



such a choice in complete detail. In other words, it still 
remains to build an actual model of dark matter within 
this framework. However, this is precisely what we shall 
do in Refs. [1,0] , and we shall verify there that our spe- 
cific models satisfy all known collider, astrophysical, and 
cosmological constraints. We thus conclude that the dy- 
namical dark-matter framework can indeed serve as a 
viable alternative to the standard paradigm of a single, 
stable, dark-matter particle, and that dynamical dark 
matter therefore has a legitimate place alongside other 
approaches to dark-matter physics. 
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Appendix: Intra-ensemble decays 

Throughout this paper, we have implicitly assumed 
that each component of our dark-matter ensemble prefer- 
entially decays directly into one or more Standard-Model 
states rather than into another, lighter component within 
the ensemble. In other words, we have been assuming 
that the decays associated with the widths Ti take the 
direct extra- ensemble form (pi SM, and that such di- 
rect decays dominate over all possible intra-ensemble de- 
cays which might produce other dark-matter components 
among their end-products. However, it is easy to gener- 
alize our overall framework to include cases in which this 
assumption is relaxed. 

Towards this end, it proves useful to start by rewriting 
Eq. ([3]) in a more useful form. Recall that f2j = Pi/pa-a, 
where pi is the energy density associated with our os- 
cillating <j>i field and where jO cr jt = 3MpH 2 , where Mp 
is the reduced Planck mass. Eq. §3§ can therefore be 
rewritten as 

Pi + iZH + Tfipi = . (A.l) 

However, because this energy density pi is entirely as- 
sociated with the coherent oscillations of the (zero- 
momentum modes of the) scalar <f>i field, it is possible 
to repackage this energy density in terms of an effective 
number density = pi/mi. We then find 

hi + (3H + Ti)ni = . (A.2) 

Indeed, Eq. (|A.2|) describes the evolution of the number 
densities associated with each of the oscillating compo- 
nents in our ensemble under the assumption that the only 
decays available for these components are direct decays 
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into Standard- Model states, i.e., <pi — > SM, with widths 

Let us now consider what happens if we introduce an 
additional set of intra-ensemble decays of the form 



(A.3) 



Here the a-index labels the specific decay channel, and 
are non-negative integers describing the multiplic- 
ities of the <fij particles produced in this decay channel 
(each of which may be assumed to have rrij < mi). Like- 
wise, collectively represents any fields outside our 
dark-matter ensemble (potentially including Standard- 
Model fields) which may also happen to be produced in 
this decay process. We shall let F^ denote the width 
associated with the decay in Eq. (|A.3|) . 

The inclusion of such additional decay channels into 
our discussion leads to two additional effects on the time- 
evolution of the number densities rij. First, there will be 
an additional decline in m due to the new decay channels 
for <f>i which are now available. However, there is also 
the possibility of an increase in rii due to the production 
of <f>i particles from the decays of presumably heavier 
components <pj within the ensemble. Indeed, we find that 
Eq. (|A.2[) is now replaced with the coupled system of 
differential equations 



Given the result in Eq. (|A.4[) . it might seem at first 
glance to be a straightforward exercise to obtain a cor- 
responding set of coupled differential equations for the 
energy densities pi and the abundances Indeed, all 
that would seem to be necessary would be to start with 
Eq. (|A.4[) in place of Eq. (IA.2p and essentially reverse the 
process that led from Eq. (|A.I[) to Eq. (|A.2[) . However, in 
going from Eq. (jA.II) to Eq. (|A.2[) we needed to assume 
that all of the energy density pi was in the form of coher- 
ent zero- momentum mode oscillations of the field (pi , and 
this will no longer be true when intra-ensemble decays are 
possible. Indeed, the daughters (pi which are produced 
through such intra-ensemble decays are literal particles 
— they have their own momenta and energies which are 
governed by the kinematics of the specific intra-ensemble 
decays which produced them. As a result, while it is 
still valid to discuss the time-evolution of a total num- 
ber density ni as in Eq. (|A.4|) . we cannot simply identify 
Pi = rami in order to obtain a corresponding set of equa- 
tions for the energy densities pi or abundances £li. 

In order to handle this calculation correctly, it is first 
necessary to express the relations in Eq. (|A.4|) in terms 
of the phase-space distributions fi{\pi\,t) [or equivalently 
fi(Ei, t), simply denoted fi] associated with each field (pi. 
This will essentially yield a Boltzmann equation which 
describes the time-evolution of these distributions. To 
do this, we observe that in general we may write 



3H 



(a) 



E £<r< a) U 



(A.4) 

In general, the solutions to Eq. (|A.4I) can exhibit a 
number of striking behaviors. Not only can there be di- 
rect decays into Standard-Model states, as before, but 
there can also be "cascade" decays that take place en- 
tirely within the dark-matter ensemble, from heavier 
states down to lighter states. Indeed, a given state can 
also decay directly into Standard-Model states at any 
point along the cascade. As a result, a particularly rich 
and subtle phenomenology can easily ensue depending 
on the relations between Ti and r^ a ' , with different por- 
tions of the ensemble exhibiting different particle-decay 
patterns in a manner reminiscent of the vacuum-decay 
patterns studied in Ref. (T^]. The corresponding values 
of rii can then alternatively rise and fall as time evolves. 

The possibility of such dark-to-dark intra-ensemble de- 
cays also allows even more striking features to emerge. 
For example, once a given heavy state (pi "turns on" , it 
can potentially decay to lighter states <pk which, because 
of their relative lightness, have not yet turned on. Thus, 
in this way, we see that a given component <pk within 
our ensemble can simultaneously contribute to dark mat- 
ter (in the form of daughter particles from the decays of 
heavier dark-matter components (pi)] to "dark radiation" 
(if the momenta of these <pk daughter particles are large 
compared to their masses); and to dark energy (in the 
form of the energy still trapped in the overdamped field 



d 3 & 

(2tt)* 



fi 



Ti = 



^ J [dw a (l ± /„)] [dn b (l ± f b )] ■ ■ ■ \M\ 2 , 



(A.5) 



where Ti is the width for a generic decay of the form 
(pi — > a + b + where M is the corresponding matrix 
element (including an implicit Dirac <5-function to en- 
force momentum conservation), where the momentum- 
integration measures are given by dir = gd 3 p/[(2ir) 3 2E] 
with g signifying the number of associated degrees of free- 
dom, and where the ± signs are chosen positive for bosons 
and negative for fermions. The left side of Eq. (|A.4|) then 
takes the form 



d 3 Pr 
(2nf 



fi+ 3ir + r i + £r^ U 



(A, 



and the right side of Eq. (|A.4[) takes the form 



j,a J I kjti 



[<M1 + fi)fi* [dn x (l±f x )]\M\ 2 (A.7) 

where Ai is the matrix element for the decay <pj — > 
12k + X^ a \ Equating the d 3 pi integrands in 
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Eqs. (IA.6I) and (IA.7I) then yields the result 
fi+^H + Ti + ^T^fi = 

[dn(l + fi)] N >?- 1 [d*x(l±fx)]\M\ 2 

(A.8) 

where we have recognized that although there are in prin- 
ciple different ways of identifying an integrand with 
respect to d 3 pi within Eq. (|A.7|) , each yields the same re- 
sult and therefore any choice will suffice. 

Eq. (|A.8[) is a set of coupled differential equations for 
the phase-space distributions ft. Indeed, while I\ is in- 
dependent of the the quantity 1^ has a hidden de- 
pendence on all fj for which N+j ^ 0. However, despite 
the complexity of this system of equations, our remain- 
ing task is conceptually easy: we simply begin with the 
distributions 

fi(Wi\,to) = 47r 3 m t (^ ) )H 3 (p l ) (A.9) 

at the time to when our abundances are initially es- 
tablished, and then use the results in Eq. ()A.8|) in or- 
der to evolve these distributions forward in time. In- 
deed, the initial distributions in Eq. (|A.9[) reflect noth- 
ing more than the assertion that the original state of 



our ensemble consists of fields 4>i whose zero-momentum 
modes are displaced by some amount (j)^ from their 
minima in field space, with a resulting energy density 
given by pi = \m\ (0-°' ) ) 2 . While the time-evolution of 
the effective number densities rij can then be obtained 
from Eq. (|A.8[) by integrating this equation with re- 
spect to the measure J d 3 pi/(2-K) 2 [thereby reproducing 
Eq. (|A.4|) ] . the time-evolution of the corresponding en- 
ergy densities pi can be obtained from Eq. (|A.8|) by in- 
tegrating Eq. (|A.8|) with respect to the alternative mea- 
sure f[d 3 p i /(2n) 3 ]E i = J[d 3 p l /{2n) 3 y\p l \ 2 + mi. This, 
then, provides us with the desired coupled differential 
equations for the energy densities pi, from which it is 
then trivial to obtain the corresponding equations for the 
abundances . 



Needless to say, there are many caveats which must 
be borne in mind when applying this formalism. First, 
in general we must require that such intra-ensemble de- 
cays not produce daughter particles with great momenta, 
for then our dark matter would not be sufficiently cold. 
Likewise, although we have included the possibility of 
intra-ensemble decays in the above analysis, we have dis- 
regarded the possible contributions from scattering pro- 
cesses which also potentially involve our ensemble compo- 
nents. Indeed, this is generally an excellent approxima- 
tion for gravitons, moduli, axions, and other fields which 
are very weakly coupled. We have also disregarded the 
effects of inverse decays. 
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